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SUMMARY 

This  report  describes  a program  for  predicting  infrared  radiation 

from  rocket  exhaust  gases  to  a point  using  a statistical  band  model  with 

' / 

a wavenumber  interval  o-T.25  cm.  Absorption  coefficient  data  for  water 
vapor,  carbon  dioxide,  carbon  monoxide,  hydrogen  chloride,  hydrogen  flouride, 
and  soot  (carbon  particles)  are  included  in  the  program  data  statements. 
Provision  has  been  made  for  the  future  addition  of  one  additional  radiating 
species.  Input  for  the  program  is  gas  property  data  (temperature,  pressure, 
and  constituent  mole  fractions)  in  either  an  axisymmetric  or  three-dimen- 
sional cylindrical  coordinate  system.  The  program  can  integrate  the  ra- 
diation from  the  gas  to  a point  to  provide  overall  radiation  as  well  as 
spectral  distribution  of  radiation  at  intervals  of725  cm"1.  In  addition, 
detailed  optical  data  can  be  provided  at  each  spectral  interval  when  only 
a single  line  of  sight  through  the  gas  is  considered. 
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INTRODUCTION 


Accurate  prediction  of  infrared  radiation  from  rocket  exhaust  plumes 

has  important  applications  in  prediction  of  heat  transfer  from  the  exhaust 

to  surrounding  structure  and  in  definition  of  the  spectral  distribution  of 

radiation  which  may  be  used  for  detecting  and  tracking  rocket  vehicles. 

The  computer  program  described  here  Is  a modification  of  one  originally 

* 

developed  by  NASA/MSFC  to  predict  radiant  transfer  to  the  base  region  of 
rocket  vehicles.  It  was  necessary  to  use  band  model  prediction  methods 
with  relatively  narrow  spectral  intervals  for  the  exhaust  plume  thermal 
radiation  to  provide  proper  representation  of  the  inhomogeneous  plumes. 

As  a result,  the  program  also  provides  spectral  data  which  may  be  useful  in 
predicting  detectability  of  rocket  powered  vehicles. 

The  primary  modification  to  the  original  NASA/MSFC  program  was  the 
addition  of  hydrogen  chloride  and  hydrogen  fluoride  radiation  calculations 
to  those  already  existing  for  water  vapor,  carbon  dioxide,  carbon  monoxide, 
and  carbon  particles.  In  addition, several  minor  changes  were  made  to  the 
program  to  correct  small  errors,  incorporate  some  more  recent  reconmenda- 
tions  of  line  half-width  parameters,  and  adapt  the  program  to  the  CDC 
6600  computer. 

A description  of  the  analytical  methods  used  and  the  general  program 
operation  is  presented  in  the  body  of  the  report  and  detailed  information 
such  as  a definition  of  terms,  input/output  formats,  a program  listing, 
and  a sample  problem  are  included  in  appendices. 

* The  program  modified  was  Job  Number  353887  obtained  from  Computation 
Laboratory,  George  C.  Marshall  Space  Flight  Center. 
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TECHNICAL  DISCUSSION 

The  computer  program  uses  Input  gas  property  data  (temperature,  pressure 
and  gas  constituent  mole  fractions)  along  with  geometric  data  locating  a 
point  of  interest  relative  to  the  gas..  Lines  of  sight  from  the  point  of 
interest  are  directed  through  the  gas  at  intervals  specified  in  the  input. 
Then  radiation  along  these  lines  of  sight  is  computed  and  the  results  are 
spatially  integrated  to  produce  the  resultant  radiant  flux  (watts/cm2)  at 
the  point  of  interest.  The  discussion  presented  in  this  section  will 
describe  the  analytical  basis  for  the  solution  and  the  numerical  methods 
used. 

DESCRIPTION  OF  THE  PROBLEM 

The  geometry  used  to  describe  the  problem  consists  of  two  coordinate 
systems  as  illustrated  in  Figure  1.  The  gas  properties  are  defined  at 
locations  in  a cylindrical  coordinate  system  relative  to  the  left-hand 
X,  Y,  Z coordinate  system.  For  axisymmetric  exhaust  plumes,  the  Z-axis 
is  the  center  of  the  cylindrical  system  for  the  gas  flow  field,  and 
properties  are  defined  at  points  located  by  r and  Z.  (Computer  programs  for 
preparing  input  tapes  for  axisymmetric  flow  fields  are  described  in 
Refs.  1 and  2).  For  three-dimensional  exhaust  plumes,  options  are  available 
to  simplify  the  flow  field  input  by  using  symmetry  which  may  exist  in  the 
flow  field.  These  options  allow  the  axis  of  the  cylindrical  coordinate 
system  to  be  moved  from  the  Z-axis  (see  Special  Option  21),  but  in  the 
general  three-dimensional,  the  gas  properties  are  defined  about  the  Z-axis 
at  points  located  by  r,  Z,  and  the  angle  n (where  n is  a positive  angle 
measured  from  the  X-axis  and  increasing  in  the  counter-clockwise  direction 
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when  viewed  from  Z = -®). 

The  point  of  interest  at  which  the  radiation  transfer  is  to  be 

computed  is  the  center  cf  a right-hand  U,  V,  W coordinate  system  In 

which  W is  the  surface  normal.  The  point  is  located  by  its  coordinates 

(X  , Y , Z ) and  the  orientation  of  the  surface  is  defined  in  general  by 
P p P ' 

specifying  the  nine  directions  cosines  defining  the  direction  of  the 
U,  V,  W axes  with  respect  to  the  X,  Y,  Z axes.  A simpler  input  method  is 
provided  in  which  only  a single  angle,  ft,  defining  the  surface  rotation  is 
specified,  but  this  method  is  limited  to  surface  elements  which  are 
oriented  as  though  they  were  on  surfaces  of  revolution  about  the  Z-axis 
(see  Special  Option  4). 

Radiant  flux  to  the  point  of  Interest  is  calculated  using  a spherical 
coordinate  system  in  which  the  inclination  of  a line  of  sight  to  the 
surface  normal  (W)  is  the  angle  6,  while  the  azimuth  angle  between  the  pro- 
jection of  the  line  of  sight  (s)  in  the  U-V  plane  and  the  U-;sxis  is  <£  (see 
Figure  1).  The  heat  transfer  equation  used  to  describe  the  radiation  to  the 
point  of  interest  in  the  spectral  region  between  wavenumbers  and  uy  is 


q/A  = 


ds  dw  d$  cos  e sin  0 de,  (1) 


where  x is  the  transmissivity  of  the  gas  and  N®  is  the  Planck  Function.  The 
equation  describes  the  radiation  from  the  gas  Integrated  over  a portion  of 
a hemisphere  (depending  on  <j>  and  e limits)  of  radius  L centered  at  the  point 
of  interest.  No  provision  is  made  for  radiation  entering  the  gas  from 
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other  sources. 

The  transmissivity  is  deterrined  from 

t (s ,(o)  = exp  X (s,w)j,  (2) 

where  X(s,u»)  is  the  optical  depth  of  the  gas  mixture.  In  the  general  case 
of  multiple  radiating  species,  the  optical  depth  Is  the  sum  of  the  optical 
depths  of  the  i radiating  species  in  the  gas  mixture 

X(s,w)  = £ X(s,w,1).  (3) 

The  method  used  for  predicting  the  optical  depth  depends  on  the 
species  considered.  Gaseous  species  are  described  by  statistical  (random) 
band  models  while  carbon  particles  are  treated  using  a model  similar  to 
Beer's  law. 

The  development  of  the  band  models  used  in  the  computer  program  have 
been  described  in  References  5 through  8 and  the  recommend  procedures  and 
data  are  summarized  in  Reference  9.  The  curve  of  growth  of  the  optical 
depth  as  a function  of  path  length,  u,  is  represented  for  collision  broadened 
(Lorentz)  lines  in  an  Isothermal  gas  by. 


X(T,p,pi,s,w,  i)  = k(T,w,i ) u(T,pi ,s)/{l+k(T,aj,1  JuCT.p^ ,s)/ 

[4YC  (T,p,Pi)/d(T,W,i)])-K 


(4) 


where  k is  the  absorption  coefficient,  yc  is  the  line  half-width,  d is  the 
line  spacing,  p is  the  pressure,  p^  Is  the  partial  pressure  of  the  radiating 
gas,  and  T is  the  gas  temperature.  The  functional  notation  used  in  Eq.  (4) 
is  cumbersome  and  will  be  eliminated  from  further  equations  except  where 
required  for  clarity.  However,  it  is  useful  at  this  point  to  clarify  the 
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definition  of  terms  as  used  in  the  analysis.  The  absorption  coefficient 
used  is  defined  for  a homogeneous  gas  by 

k = lim  [(-  In  t)/u],  (5) 

u— *o  L J 

where  u is  the  equivalent  optical  path  at  standard  temperature  and  pressure, 

u = c (pm)  [ P.(aM  273K*1  . 

L 1 atm  T(K)  J (6) 

The  absorption  coefficient  used  in  the  computer  prof ram  is  an  average 

value  representing  the  equivalent  absorption  in  a 25  cm"1  wavenumber  interval, 

so  fine  spectral  details  cannot  be  reproduced  using  the  program.  This 

spectral  averaging  over  a small  wavenumber  Interval  is  implied  In  all  the 

equations  describing  the  band  models  although  no  special  notation  Is  used  to 

denote  an  average  value  except  In  the  case  of  spatial  averaging  over  an 

inhomogeneous  path.  The  only  monochromatic  value  used  in  the  program  is 

the  Plank  function,  N°. 

w 

The  line  half-width  is  treated  as  a function  of  temperature,  partial 
pressure  of  the  radiating  constituent,  and  partial  pressures  of  all  the 
other  constituents.  Dependence  of  the  fine  structure  parameter 

= y /d  (7) 

c c 

on  wavenumber  is  accounted  for  by  the  line  spacing  d which  is  a function 
of  temperature  and  wavenumber  for  each  constituent. 

The  curve  of  growth  represented  by  Eq. (4)  has  the  limits 

X = ku  (8) 
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for  small  values  of  u (ku«4Y/d),  and 

X = (4  ku  Yc/d)  ^ (9) 

for  large  values  of  u (ku»4Y/d).  These  are  referred  to  as  the  weak 
line  and  strong  line  limits  which  are  characteristic  of  isolated  Lorentz 
lines  as  well  as  several  statistical  band  models.  The  behavior  of  Eq.  (4) 
between  the  limits  accurately  represents  a band  model  composed  of  randomly 
located  lines  in  which  the  probability  of  occurance  of  a given  line  strength 
Is  proportional  to  exp(-  line  strength).  This  Is  the  so-called  "Goody'* 
model,  and  widespread  use  of  this  model  results  not  so  much  from  any 
rational  physical  basis  as  from  the  fact  that  the  expression  for  the  curve 
of  growth  has  a simple  algebraic  form.  It  Is  also  a reasonable  approximation 
of  other  statistical  band  models.  For  example,  Eq.(4)  represents  the  curve 
of  growth  for  a spectrum  of  equal  strength  randomly  located  lines  with  a 
maximum  error  of  8 percent. 

Caution  should  be  exercised  In  using  the  statistical  band  model  repre- 
sentation to  describe  gases  with  relatively  regular  line  spacing  such  as 
HC1  and  HF.  If  the  results  obtained  do  not  appear  to  compare  well  with 
measured  data,  consideration  should  be  given  to  using  the  Elsasser  model 
which  assumes  a uniform  spacing  of  equally  strong  Lorentz  lines. 

The  optical  depth  due  to  Doppler  broadening  is  also  modeled  using 
a random  band  model.  The  expression  used  Ms  a mathematical  approximation 
(within  10  percent)  of  the  curve  of  growth  which  results  from  a random 
distribution  of  equal  strength  Doppler  lines.  For  conditions  in  which 
both  collision  and  Doppler  broadening  are  important,  a form  is  provided1 2 

1.  See  Table  1,  Eq.  (17a) 

2.  See  Table  1,  Eq.  (14a) 
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approximating  a line  having  combined  Lorentz  and  Doppler  broadening. 

In  the  discussion  of  band  models  up  to  this  point,  only  homogenous 
gases  have  been  considered.  The  models  are  applied  to  inhomogenous  gases 
through  a modification  of  the  Curtis-Godson  approximation  which  was 
originally  proposed  for  a single  line.  The  method  used  in  the  program 
replaces  ku  with  an  integrated  optical  depth  for  the  weak  line  limit 


and  averages  the  fine  structure  parameter  weighed  by  the  product  ku. 


This  method  of  averaging  has  been  shown  experimentally  to  give  good 
results  over  moderate  ranges  of  temperature  and  absorption  (Ref.  10)  , but 
it  may  be  expected  to  produce  errors  in  some  extreme  cases.  In  particular, 
when  thpre  are  large  temperature  variations  along  a line  of  sight  so  that 
radiation  from  a hot  gas  layer  must  pass  through  a relatively  massive 
cold  gas  layer,  errors  are  likely  to  be  significant.  Fortunately,  the 
isentropic  expansion  in  rocket  exhaust  plumes  produces  extremely  low 
pressures  in  low  temperature  regions,  so  the  mass  of  cool  layers  tends 
to  be  small. 

In  order  to  provide  a more  exact  representation  of  the  effects  en- 
countered with  drastic  temperature  variations  as  well  as  a simpler  method 
for  less  demanding  conditions,  two  forms  of  band  models  are  used  in 
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the  computer  program.  The  more  exact  representation  uses  a superposition 
of  a number  of  line  groups  to  account  for  the  additional  line  positions  as 
higher  vibrational  states  are  excited  at  high  temperatures,  while  the  simpler 
representation  neglects  the  appearance  of  new  lines.  The  existing  program 
nomenclature  of  Model  3 and  Model  3A  for  these  two  methods  will  be  re- 
tained in  this  discussion  without  attempting  to  describe  the  history  of  the 
numbering  system.  All  the  relations  comprising  the  two  models  are  listed 
in  Table  1,  and  the  required  line  width  parameters  are  presented  in  Table  2. 

The  9 term  in  Eq.  (20)  (Table  1)  represents  the  energy  associated  with 
a vibrational  transition,  so  the  exp(-0n/T)  term  represents  a measure  of  the 
probability  of  shifted  line  positions  (due  to  vibrational  transitions)  as  a 
function  of  the  gas  temperature.  The  form  of  the  absorption  coefficient  for 
Model  3 reduces  to  that  for  Model  3A  as  0 goes  to  infinity,  and  the  1/d 
values  for  the  two  models  are  formulated  to  give  the  same  optical  depth  in 
the  strong  line  limit.  If  T is  somewhat  less  than  0,  the  usual  result  Is 
that  radiance  predicted  by  Model  3 is  less  than  Model  3A.  This  is  caused  by 
the  reduction  in  optical  depth  due  to  the  reduced  line  density  of  Model  3 
being  more  significant  than  the  increases  due  to  the  added  line  groups. 

Although  Table  1 indicates  that  Model  3 consists  of  an  infinite 
number  of  line  groups,  only  a few  groups  are  required  in  practice.  Inspection 
of  Eq.  (20)  (Table  1)  indicates  that  the  exponential  term  will  become 
negligible  compared  to  the  first  term  at  small  values  of  n if  the  gas 
temperature  is  within  the  range  of  the  absorption  coefficient  data  (T<3000K). 
Because  of  this,  the  computer  program  has  been  arranged  to  handle  only  10 
line  groups  (n  = 0 to  9)  to  reduce  storage  requirements. 

The  functions  used  to  relate  1/d  and  ]/dQ  are  derived  by  equating  the 
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TABLE  2 

VALUES  FOR  THE  COLLISION 
BROADENING  LINE  WIDTH  PARAMETERS  (a) 


Molecule 

0) 

Broadener 

(J) 

(yi ,j^273 
cm"  atm"1 

Q 

^Yi.i*273 
cnr  atnf 1 

ni.i 

H20 

H20 

0.44 

1.0 

N2 

0.5 

02 

0.04 

0.5 

Ha 

(0.05) 

0.5 

C02 

0.12 

0.5 

CO 

(0.10) 

0.5 

Other  (b) 

(0.05) 

0.5 

co2 

C02 

0.09 

rm 

0.01 

1.0 

H20 

(0.07) 

H 

Ha 

0.07 

0.5 

o. 

0.055 

KB 

H2 

0.08 

EH 

CO 

(0.06) 

EH 

Other  (b) 

(0.08) 

■ 

CO 

CO 

0.06 

MH9 

0.0 

1.0 

H20 

(0.06) 

1 

C02 

(0.07) 

1 

Ha 

0.06 

1 

Ha 

0.06 

'M 

02 

0.05 

0.5 

Other  (b) 

(0.06) 

0.5 

HC1 

HC1 

mm 

■ 

0.15 

1.0 

Other  (b) 

111 

19 

HF 

HF 

mm 

19 

0.45 

1.0 

Other  (b) 

Ia9 

19 

NOTE:  (a)  Values  in  parenthesis  are  estimated. 

(b)  Values  for  other  gases  listed  here  are  used  in  the  program,  but  they 
may  be  modified  using  input  Card  25. 
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optical  depth  of  Model  3 and  Model  3A  In  the  strong  line  limit, 
gi  ves 


This 


(22) 


which,  with  suitable  manipulation,  becomes 


1/d  = l/dn  | [l+exp(-6/2T)]y/  [l -exp (-8/2T)]  | 


(23) 


for  fn  = 9n  = 1*  for  the  case  °f  fn  = n+l  an£l  9n  = 1»  Eq.  (22) 
may  be  approximated  within  3 percent  by 


1/d  - l/do 


(24) 


The  data  used  in  the  computer  program  for  the  band  model  parameters 
k,  Y,  and  1/d  are  based  on  a combination  of  experiments,  analyses,  and 
estimates  based  on  the  expected  characteristics  of  the  gaseous  radiators. 

The  collision  half-width  parameters  listed  in  Table  2 are  a combination  of 
experimental  and  estimated  values  (Ref.  9).  Values  of  k and  1/d  for  water 
vapor  in  the  areas  of  primary  importance  were  determined  by  fitting  the 
band  model  for  collision  broadening  to  experimental  measurements  (Refs.  7 
and  11).  In  regions  where  no  data  were  obtained,  estimates  of  the  radiation 
behavior  were  made.  Values  of  k and  1/d  for  the  remaining  gases  are  based 
primarily  on  theoretical  predictions  of  the  detailed  line  structure  which 
is  fitted  by  the  band  model  (Refs.  5,  6,  and  9).  (It  should  be  noted  that 
temperature  and  wavenumber  dependent  values  for  1/d  for  CO  are  available 
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(Ref.  9),  but  these  have  not  been  Incorporated  In  the  computer  program). 

Investigations  of  the  carbon  particles  occuring  in  rocket  exhausts 
(Refs.  3 and  4)  indicated  that  the  particles  were  so  small  that  scattering 
effects  were  not  significant  for  problems  of  heat  transfer  from  exhaust 
plumes.  Because  of  this,  scattering  has  been  neglected  in  the  computer 
program,  and  carbon  particles  are  treated  as  continuum  radiators  following 
a Beer's  law  model  using  experimentally  determined  values  of  the  absorption 
coefficient.  The  optical  depth  for  carbon  particles  is 

s 

v/  kcpcds‘ 

*o 

where  k is  the  absorption  coefficient  per  unit  density  and  path  length 

(cn2/gm),  p is  the  local  carbon  particle  density  per  unit  gas  volume 
0 

(gm/cm3),  and  s is  the  path  length  (cm). 

Experimental  values  of  the  absorption  coefficient  (illustrated  in 
Fig.  2)  could  only  be  determined  over  the  wavelength  range  from  1 pm  to 
4 pm  (2500  cm  1 to  10,000  cm"1  wavenumbers).  Since  serious  errors  could 
result  in  heat  transfer  by  limiting  calculations  to  this  range,  the 
data  were  extrapolated  following  the  trends  in  the  data  to  cover  a range 
from  1 pm  to  10  pm  (1000  cm"1  to  10,000  cm"1).  The  absorption  coefficient 
data  is  stored  in  the  program  as  polynomials  which  are  functions  of 
wavenumber  for  each  the  temperatures  indicated  in  Fig.  2. 
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Figure  2.  Results  of  the  experimental  determination  of  the  carbon 
absorption  coefficients.  (Ref.  3) 
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NUMERICAL  METHODS  OF  SOLUTION 

The  numerical  techniques  used  in  the  program  are  relatively  straight- 
forward. The  gas  properties  are  determined  at  points  along  lines  of  sight 
by  interpolation  in  the  input  gas  data,  and  the  radiation  along  each  line  of 
sight  is  computed  using  summations  to  represent  the  integration  procedure 
described  in  the  previous  section. 

Computations  along  the  lines  of  sight  begin  with  the  line  represented 
by  the  lower  limits  of  the  0 and  ^integration  (6  * 0j  + A8/2  and  4>  = <f>^  + 
Acf>/2) . The  sequence  of  processing  the  lines  carries  through  all  values  of 
<j>  before  starting  with  the  next  interval  of  0 at  the  minimum  value  of  <j>. 

In  calculating  the  positions  on  each  line  of  sight,  the  direction  cosines 
of  the  line  with  respect  to  the  X,  Y,  and  Z axes  are  determined  first,  then 
the  intersection  of  the  line  of  sight  and  the  blocking  circle  positions 
are  checked  to  determine  if  the  length  of  the  line  of  sight  needs  to  be  ad- 
justed due  to  occlusion  by  one  of  the  circles.  Next,  increments  AX,  AY,  and 
AZ  corresponding  to  the  input  value  of  As  are  calculated,  and  the  coordinates 
of  the  center  of  the  first  element  on  the  line  of  sight  is  completed  (X  = Xp  + 
AX/2,  Y = Yp  i AY/2,  and  Z = Zp  + AZ/2).  Coordinates  of  each  subsequent 
point  are  computed  by  successively  incrementing  the  coordinates  by  AX,  AY, 
and  AZ. 

As  the  coordinate  of  each  point  is  computed,  it  is  tested  to  determine 
if  it  is  within  the  flow  field.  If  it  is  not,  the  temperature  at  that 
point  is  set  equal  to  zero  and  the  next  point  is  considered.  If  it  is  in 
the  flow  field,  a modified  form  of  linear  interpolation  is  used  to  determine 
the  gas  properties.  The  procedure  used  for  axisymmetric  flow  fields  is 


19 


REMTECH  INCORPORATED 


RTR  013-1 


illustrated  below  (subscript  w is  the  point  at  which  any  property  P is 
to  be  computed) 


ZR  = (Z  w - Z1)/(Z2  - Z,) 

Rwmax  = Rlmax  + ZR^R2max  " Rlmax^ 
R1  ~ ^max^w/^wmax^ 

R2  _ R2max^Rv/Rwmax 
pw  = P1  + zp'(p2  “ pl> 


Z 


If  the  shock  position  is  not  specified,  the  properties  P-j  and  P,,  are  those 
computed  at  R^  and  R2  using  linear  interpolation  between  the  points  on  either 
side  of  R.|  and  R^.  If  the  shock  position  is  specified  (see  Special  Option  10), 
the  interpolation  procedure  is  modified.  For  the  case  illustrated  in  which 
the  shock  position  is  outside  the  point  under  consideration,  only  properties 
at  points  inside  the  shock  would  be  used  for  interpolation.  Therefore,  the 
properties  at  points  R-]a  and  R^  (the  last  two  points  inside  the  shock  lo- 
cation) are  linearly  extrapolated  to  the  radius  R^ . This  gives 

P1  * PU  +<P1I>  - Pla)(Rl  - Rla>/<P'lb  - Rl>> 

Similarly,  if  Rw  were  just  outside  the  shock,  R-j  would  be  treated  normally, 
but  the  properties  of  R2  would  be  extrapolated  from  R2a  and  Rg^  (the  first 
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two  points  following  the  last  point  inside  the  shock). 

In  the  three-dimensional  flow  field  interpolation,  the  interpolation 
on  R is  carried  out  as  described  above  for  the  two  values  of  n surrounding 
the  point  to  be  defined.  Then  an  interpolation  in  the  r\  direction  is  carried 
out  in  the  same  fashion  as  that  just  described  for  R,  so  that  the  interpolated 
points  on  each  n plane  have  the  same  radius  ratio  relative  to  the  assumed  gas 
boundary.  After  properties  are  determined  at  each  point  on  each  line  of 
sight,  the  radiation  calculation  is  startad. 

The  use  of  stepwise  summations  to  represent  the  integrals  in  the  radia- 
tion calculations  saves  computer  time  and  storage  and,  when  carefully  handled, 
is  expected  to  give  accuracy  consistent  with  the  definition  of  the  gas  proper- 
ties and  band  model  parameters.  Caution  must  be  observed  in  selection  of  the 
geometric  intervals  of  integration,  so  that  results  represent  the  variations 
in  gas  properties.  The  increment  along  the  line  of  sight.  As,  should  be  chosen 
so  that  large  changes  in  transmissivity  do  not  occur  in  any  single  increment. 
Achievement  of  this  goal  can  be  judged  to  some  extent  by  evaluating  the  changes 
in  integrated  transmissivity  along  a line  of  sight  (using  NLINE  = 2 or  3).  If 
the  transmissivity  changes  greatly  between  Integration  steps  in  a region  of 
rapidly  changing  properties,  then  a smaller  step  size  should  be  tried.  A 
temperature  criteria  (Special  Option  19)  has  been  provided  to  allow  the  use  of 
a small  step  size  to  define  regions  of  rapid  change  without  making  unnecessary 
calculations  in  regions  of  relatively  constant  properties. 

An  evaluation  of  the  numerical  accuracy  of  the  angular  integration  can 
be  made  by  comparing  the  output  exact  value  of  the  total  "shape  factor"  with 
the  numerical  value  (where  "shape  factor"  = /sin  6 cos  e d9d<f>) . However, 
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determining  if  the  angular  increments  are  small  enough  to  define  the  property 
variations  in  the  gas  has  to  be  a matter  of  judgement. 
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FUNCTIONAL  STRUCTURE 

The  program  is  divided  Into  two  overlay  subroutines  which  call  additional 
overlay  subroutines.  The  two  major  subroutines,  FLOW  and  RAD,  separate  the 
functions  of  interpolation  in  the  gas  property  data  and  computation  of  the 
radiation,  so  the  flow  field  property  data  and  the  radiation  coefficient  data 
do  not  have  to  be  stored  simultaneously.  Functional  arrangement  of  the  program 
and  the  logic  used  in  calling  the  subroutines  is  illustrated  in  Figure  3,  and 
a brief  description  of  the  function  of  each  subroutine  is  given  below. 

MAIN 

This  subroutine  reads  the  initial  control  cards  and  controls  calling  the  two 
major  subroutines. 

FLOW 

This  controls  the  subroutines  which  output  the  case  data  and  prepare  the  flow 
field  gas  property  data  for  the  radiation  calculation. 

FLOWIN 

This  subroutine  reads  the  initial  identifying  portions  of  the  flow  field  input 
and  the  input  data  necessary  for  the  flow  field  interpolation,  and  writes  the 
input  data  on  tape  10.  It  also  identifies  the  constituents  in  the  flow  field 
and  indexes  them  for  later  manipulation. 

FLOUT 

This  subroutine  reads  tape  10  to  prepare  the  initial  printed  output  describing 
the  problem  limits  and  leaves  tape  10  in  position  to  read  or  write  the  first 
set  of  interpolated  data  along  a line  of  sight. 
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FL0W2D 

? This  subroutine  prepares  gas  property  data  along  lines  of  sight  by  interpo- 

lation in  an  axi symmetric  flow  field  and  writes  these  data  on  tape  10.  When 
I the  program  returns  from  this  subroutine,  the  line  of  sight  data  on  tape  10 

I is  complete  for  the  problem  limits  specified.  This  subroutine  also  provides 

printed  output  of  the  flow  field  properties  and  prepares  a flow  field  on 
I tape  8 if  IFL0TP>0. 


FL0W3D 

This  subroutine  prepares  gas  property  data  along  lines  of  sight  by  interpolation 
in  a three-dimensional  flow  field  and  writes  these  data  on  tape  10.  When  the 
program  returns  from  this  subroutine,  the  line  of  sight  data  on  tape  10  is 
complete  for  the  problem  limits  specified.  This  subroutine  also  provides  printed 
output  of  the  flow  field  properties,  and  prepares  a flow  field  on  tape  8 if 
IFL0TP>0. 


SKIP 

The  subroutine  SKIP  renews  the  output  page  and  line  count.  It  is  called  from 
FLOUT  only. 


RAD 

This  controls  the  operation  of  subroutines  for  the  radiation  calculation  and 
output. 

AC DATA 

This  subroutine  provides  for  modification  of  the  absorption  coefficient  data 
as  required,  and  provides  printed  output  of  the  absorption  coefficient  data 
used. 
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M0DL3A 

This  subroutine  calculates  radiation  using  the  method  described  as  Model  3A. 
The  line  of  sight  data  on  tape  10  is  used  as  input,  and  printed  output  of 
the  spectrally  integrated  radiation  data  is  provided. 

M0DEL3 

This  subroutine  calculates  radiation,  using  Model  3 equations  with  the  number 
of  line  groups  n + 1 = IRAD.  It  uses  tape  10  as  input  and  provides  printed 
output  of  the  spectrally  integrated  radiation  results. 

SPCOUT 

This  subroutine  provides  for  printing  and  plotting  output  for  the  spectral 
distribution  of  the  radiation  data. 
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Figure  3.  Program  Arrangement 
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INPUT/OUTPUT 

The  program  can  be  operated  with  card  input  and  printed  output,  but 
options  using  tapes  and  providing  plotted  output  are  available  fo.r  convenience 
Band  model  data  is  included  in  the  program,  but  the  1/d  values  for  CO^  must 
be  modified  using  card  input  for  making  calculations  with  Model  3 (refer  to 
Description  of  the  Problem). 

Geometric  description  of  the  gas  property  data  can  be  input  using  a tape 
generated  by  an  axisymmetric  method-of-characteristics  program  (Refs.  1 and  2) 
The  radiation  program  can  also  generate  a flow  field  tape  for  later  use  from 
either  card  or  tape  input  with  modifications  (such  as  scaling  or  carbon 
addition)  dictated  by  the  radiation  program  input. 

The  tape  with  the  gas  property  data  ordered  along  each  line  of  sight, 
which  is  generated  normally  as  a working  tape,  can  be  saved  and  used  as  in- 
put. This  is  convenient  for  restarting  an  incomplete  calculation  or  for 
running  a problem  again  with  scaled  properties  or  reduced  integration  limits. 

Plotted  output  is  available  for  the  spectral  distribution  of  the  Integra- 
ted radiation  and  for  the  transmissivity  and  optical  properties  on  the  last 
line  of  sight  processed. 

Detailed  descriptions  of  the  input  and  output  formats  are  given  in 
Appendix  B. 
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SPECIAL  OPTIONS 

This  program  has  many  special  options.  They  will  be  briefly  described 
in  the  order  in  which  they  appear  in  the  input  data. 

1.  The  gas  property  data  (flow  field)  may  be  input  from  cards  or 
tape.  If  desired,  a previously  generated  line-of-sight  tape  may 
be  rerun  (IFL0W<1)  with  changes  in  only  output  form  or  scaling 
(KASE  =0).  However,  if  problem  limits  are  to  be  reduced,  KASE 
must  match  the  case  number  used  when  generating  on  tape  10. 

2.  An  input  flow  field  may  be  modified  by  the  addition  of  a specified 
carbon  mole  fraction  (see  Op i ton  6)  and  then  output  on  another  tape 
which  may  or  may  not  contain  additional  cases.  This  allows  flow 
fields  to  be  put  on  tape  from  cards,  modified  and  preserved  to 
simplify  future  runs,  or  consolidated  to  reduce  the  number  of  tape 
being  used. 

3.  The  radiation  calculation  may  be  performed  from  a line-of-sight 
tape  previously  generated,  or  from  either  axisymmetric  or  three- 
dimensional  gas  property  input  data. 

4.  The  plane  of  the  point  of  interest  in  the  radiation  calculation 

can  be  defined  either  by  nine  direction  cosines  or  by  a single  angle 
of  inclination.  When  the  angle  of  inclination  (OMEGA)  is  specified, 
the  Z-axis  must  be  in  the  U-W  plane.  OMEGA  is  the  angle  measured 
from  the  +Z-axis  t--  the  +W-axis.  With  OMEGA  equal  0,  the  positive  U 
axis  intersects  and  is  perpendicular  to  the  Z-axis.  Positive  OMEGA 
ratation  is  clockwise  about  the  V-axis  when  looking  in  the  positive 
V direction. 

5.  Up  to  50  blocking  circles  may  be  used  to  simulate  structure  which 
may  block  the  radiation  lines  of  sight.  The  blocking  circles  may  be 
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either  "disks"  or  "holes",  and  each  circle  Is  located  by  the  coordinates 
of  Its  center  and  Its  radius.  When  a line-of-sight  passes  through  a 
"disk"  or  outside  a "hole"  it  is  terminated  at  the  blocking  circle. 

The  blocking  circles  are  restricted  to  being  parallel  to  the  X-Y 
plane. 

6.  Carbon  particles  (soot)  may  be  added  in  three  ways  to  flow  fields 
which  do  not  contain  them.  The  carbon  mole  fraction  (treating  it 
as  a gas  with  mole  weight  of  12.01)  may  be  added  as  (1)  a constant, 

(2)  a function  of  radial  position,  or  (3)  as  a function  of  radial 
and  axial  position. 

7.  The  radiation  calculation  results  may  be  output  in  any  of  four  forms 
with  varying  degrees  of  detail.  These  methods  are  explained  in  the 
input  instructions. 

8.  The  program  is  arranged  to  display  several  forms  of  spectral  output 
including  plots.  See  input  instructions  for  further  definition. 

9.  To  reduce  output  bulk,  abbreviated  forms  of  both  flow  field  and 
absorption  coefficient  output  are  provided.  In  the  abbreviated 
form,  only  the  first  and  the  last  entry  in  each  block  of  data  is 
printed  to  assure  proper  input  operation. 

10.  In  order  to  define  the  discontinuity  in  gas  properties  caused  by 
a shock  wave,  an  option  is  provided  to  input  the  exact  shock  posi- 
tion. This  allows  the  flow  field  interpolation  procedure  to  identify 
the  position  of  the  point  being  defined  relative  to  the  shock  wave,  so 
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that  the  interpolation  may  be  modified  when  required.  The  result  is 
a sharp  break  at  the  shock  position  between  hot  and  cold  gas  rather 
than  a smeared  out  warm  region.  This  option  may  use  a shock  location 
defined  in  the  flow  field  (ISHOCK  = -1)  or  input  separately  (ISHOCK  = 1) 
If  this  option  is  to  be  used  there  must  be  at  least  two  radial  points 
on  each  side  of  the  shock. 

11.  The  program  does  all  calculations  in  the  MKS  system  of  units,  but  the 
geometric  and  gas  property  data  may  be  input  in  either  the  English 
(inches,  R,  lbf/ft2)  or  the  MKS  (cm,  K,  atm)  system. 

12.  The  flow  field  dimensions  and  properties  may  be  scaled,  in  order  to 
perform  parametric  studies  or  convert  dimensions  to  the  proper  units. 
When  using  previously  generated  line  of  sight  tape  (IFL0W<1 ),  the  scale 
factors  will  be  multiplied  by  scale  factors  used  in  generating  the  tape. 

13.  If  it  is  desired  to  load  less  than  the  entire  available  input  flow  field 
the  position  of  the  initial  and  final  Z-planes  may  be  specified. 

14.  The  radiation  calculation  may  be  performed  using  either  Model  3 or 
Model  3A  radiation  equations  (see  technical  discussion). 

15.  The  program  contains  the  absorption  coefficient  data  which  will  normally 
be  required  for  Model  3A,  but  provisions  are  made  for  modifying  this 
data  without  changing  it  permanently  if  It  is  necessary  for  a particular 
problem.  The  data  contained  in  this  program  consists  of  1620  entries 
for  constituent,  parameter,  or  wavenumber  with  7 temperature  values  for 
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each  entry.  The  arrangement  is: 


Constituent 

Parameter 

Position 
in  Table 

Lower  Value 

of  V 

Upper  Value 

of  V 

Av 

h2° 

k 

1-439 

50 

11000 

25 

1/d 

440-878 

50 

11000 

25 

CO, 

k 

879-1009 

500 

3750 

25 

1/d 

1010-1140 

500 

3750 

25 

CO 

k 

1141-1194 

1025 

2350 

25 

HC1 

k 

1195-1284 

1000 

3225 

25 

1/d 

1285-1374 

1000 

3225 

25 

HF 

k 

1375-1497 

1400 

4450 

25 

1/d 

1498-1620 

1400 

4450 

25 

The  analytical  expression  listed  for  H20  1/d  in  Table  1 can  be  used 
rather  than  the  tabulated  values  by  setting  ICOEFF  =-l.  (In  normal 
usage  this  will  not  usually  be  required,  and  the  logic  which  allows  the 
analytical  expression  in  Model  3 and  Model  3A  could  be  removed  to  save 
time.  On  the  other  hand,  if  storage  becomes  a problem,  the  tabulated 
values  could  be  eliminated  in  favor  of  the  analytical  expression.) 

When  modifying  the  tabulated  values  of  k or  1/d  (i.e.,  to  input  l/dQ 
used  in  Model  3 for  C02  or  to  change  the  wavenumber  range  or  inter- 
val), two  restrictions  must  be  observed: 

(a)  The  total  number  of  entries  must  not  exceed  1620  and  the  number 
of  temperature  values  for  each  entry  must  not  exceed  7. 

(b)  If  both  k and  1/d  are  tabulated  for  a constituent,  the  wave- 

number  range  for  both  parameters  must  be  identical.  For  example 

if  1/d  is  modified  for  CO  , the  input  tabulated  data  must  contain 

2 

132  entries  starting  at  w=  500  with  Aof25.  However,  if  both  k 
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and  1/d  are  modified,  a wavenumber  range  different  from  that  existing 
in  the  program  may  be  used. 

16.  The  absorpiton  coefficients  may  be  scaled  in  order  to  perform  parametric 
studies  or  convert  dimensions  to  the  proper  units  (cnT^yp* 

17.  In  order  to  provide  more  uniform  spatial  integration,  the  parameter  TANGLE  = 
cos  9 sin  0 A0Acf>  may  be  specified  (card  4).  The  value  of  TANGLE  will  be 
used  to  compute  a value  of  DPHI  to  be  used  for  the  Integration  increment 

at  each  value  of  0.  However,  if  it  is  smaller  than  the  input  DPHI,  the 
input  DPHI  will  be  used  as  the  integration  increment. 

18.  A line  of  sight  tape  may  be  prepared  by  the  program  without  making  a 
radiation  calculation  by  setting  DENU  = 0. 

19.  A temperature  step  criteria,  TDIFF,  is  provided  to  vary  the  geometric 
Integration  step,  As,  depending  upon  the  variations  of  temperature  in 
the  flow  field.  With  TDIFF  >0,  the  flow  field  is  interpolated  using 
steps  of  As  and  summing  the  property  values  until  the  temperature  change 
from  the  first  point  in  the  integration  increment  is  approximately  equal 
to  TDIFF.  When  this  occurs,  the  property  values  are  averaged  and  a 
step  in  the  radiation  calculation  is  computed. 

20.  An  option  is  provided  to  exclude  insignificant  absorption  coefficients 
from  the  radiation  calculation  to  save  time.  If  an  absorption 
coefficient  is  less  than  or  equal  to  ACMIN,  it  will  be  neglected. 

21.  The  form  of  three-dimensional  gas  property  data  has  been  made  as 
flexible  as  possible  to  reduce  the  number  of  data  points  by  taking 
advantage  of  symmetry  such  as  that  resulting  from  a uniform  arrange- 
ment of  multiple  rocket  engines.  The  form  used  is  presented  to 
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the  program  by  the  value  of  three  input  parameters:  HANG,  CTCD, 
and  ISYMSC.  Several  typical  examples  are  presented  in  Fig.  4 to 
clarify  the  form  of  the  gas  property  coordinate  system  and  the  use 
of  the  three  input  parameters.  In  the  general  case  (Fig.  4(a)),  the 
angle  about  the  Z-axis  over  which  property  data  is  provided  (HANG)  is 
a full  360  degrees.  With  symmetry  about  the  XZ  plane,  HANG  is  reduced 
to  180°  (Fig.  4(b)),  and  addition  of  symmetry  about  the  YZ  plane  allows 
a reduction  of  HANG  (and  n)  to  90  degrees.  If  it  Is  desired  to  offset 
the  axis  of  the  cylindrical  coordinate  system  for  the  gas  properties 
as  in  Fig.  4 (d)  , the  offset  distance  is  CTCD.  In  using  this  option, 
the  axis  of  the  cylindrical  coordinate  system  must  be  in  the  XZ  plane 
parallel  to  the  Z-axis,  and  r\  and  R must  be  chosen  so  that  the 
gas  properties  are  bounded  by  the  XZ  and  YZ  planes.  This  technique 
can  be  extended  to  smaller  symmetric  sectors  as  in  Fig.  4(e).  For 
cases  in  which  a center  engine  is  used,  it  may  be  desirable  (although 
not  technically  correct)  to  assume  symmetry  about  a plane  between  the 
center  and  outboard  engine  as  in  Fig.  4(f).  In  this  case  ISYMSC  is 
set  equal  0,  and  n and  R must  be  chosen  so  that  the  gas  properties 
are  bounded  by  CD  as  well  as  the  radial  planes  of  synmetry. 

In  all  cases  illustrated  in  Fig.  4,  the  maximum  range  of  n is 
Irdi cated.  However,  if  the  gas  terminates  before  the  upper  limit 
of  n,  it  is  not  necessary  to  carry  the  flow  field  description  to 
the  upper  limit. 
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Axis  of  Cylindrical  Coordinates  Parallel  to  Z 


Rocket  Nozzle 
Exit  Plane-7 


HANG 

Comments 

360° 

0 

o 

oo 

90° 

90° 

R bounded  by 
YZ  plane 

45° 

R bounded  by 
AZ 

45° 

R bounded  by 
CD  and  AZ 

(a)  General 

(b)  Symmetric  about  XZ  plane 

(c)  Symmetric  about  XZ  and  YZ 
planes 

(d)  Symmetric  about  XZ  and  YZ 
planes  and  flow  field  off 
center 

(e)  Symmetric  about  XZ,  YZ,  AZ, 
and  BZ  planes 

(f)  Symmetric  sectors  as  In  (e) 
plus  symmetry  of  regions  E 
and  F about  CD 


Fig.  4.  Examples  of  Techniques  for  specifying  three  dimensional  flow  fields. 
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RESTRICTIONS  AND  LIMITATIONS 

This  program  has  many  liritations  caused  by  storage,  but  methods  have 
been  incorporated  which  remove  some  of  the  limits  imposed  on  previous  gaseous 
radiation  programs. 

The  general  limitations  on  the  flow  field  subroutines  are  as  follows: 
Maximum  number  of  constituents  to  be  stored  is  6 
Maximum  number  of  Z-planes  used  for  the  addition  of  carbon  (ICARB>0) 
is  10 

Maximum  number  of  rj-values  used  for  the  addition  of  carbon  (ICARB>0) 
per  Z-plane  Is  7 

Maximum  number  of  radial  points  used  for  the  addition  of  carbon 
(ICARB>0)  is  10 

Maximum  number  of  Z-planes  at  which  shock  radii  may  be  specified  when 
shock  radii  are  added  (ISHOCK  = 1)  is  20 

Maximum  number  of  n-values  at  each  Z-plane  for  which  shock  radii  may 
be  specified  is  13 

Maximum  number  of  blocking  circles  is  50 
In  addition  to  these  limitations,  each  flow  field  interpolation  subroutine 
(FL0W2D  and  FL0W3D)  has  specific  limitations.  In  each  of  these  sub- 
routines, the  number  of  Z-planes  which  may  be  used  is  unlimited,  since  the 
program  is  arranged  to  handle  a portion  of  the  flow  field  at  a time. 

In  FLQW2D,  the  maximum  number  of  radial  points  permitted  is  100.  If 
the  number  of  points  exceeds  100,  points  wl1!  be  skipped  to  reduce  the  data 
to  less  than  100  points.  The  maximum  number  of  Z-planes  which  is  permitted 
in  storage  is  10,  but  the  program  will  automatically  subdivide  the  flow 
field  so  that  any  number  of  Z-planes  may  be  used.  However,  the  length  of  the 
flow  field  which  can  be  used  is  regulated  indirectly  by  a limit  of  300  points 
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along  any  line  of  sight.  If  this  limit  is  exceeded,  a diagnostic  message  is 
printed,  and  the  calculation  is  stopped. 

In  FL0W3D,  the  maximum  number  of  radial  points  permitted  is  30.  If 
the  number  of  points  exceeds  30,  points  will  be  skipped  to  reduce  the  data 

i 

to  less  than  30  Doints.  The  maximum  number  of  n-values  at  any  Z-plane  is 
13.  If  this  number  is  exceeded,  the  program  will  stop  and  print  a diagnostic 
message.  Only  3 Z-planes  are  permitted  in  storage  at  any  time,  but  the 
program  automatically  subdivides  the  flow  field  so  that  any  number  of 
Z-planes  may  be  used.  The  maximum  number  of  points  permitted  along  any  line 
of  sight  is  300.  If  this  limit  is  exceeded, a diagnostic  message  is  printed, 
and  the  program  stops.  Geometry  restrictions  are  described  in  Special  Option  21. 

In  the  radiation  calculation  subroutines,  the  limits  on  the  number  of 
constituents  and  number  of  points  allowed  on  a line  of  sight  are  compatible 
with  the  flow  field  subroutines.  The  primary  limitations  of  interest  in 
the  radiation  subroutines  are  those  pertaining  to  storage  of  spectral  and 
constituent  radiation  data. 

Space  is  provided  for  a table  of  coefficients  with  up  to  1620  entries 
and  up  to  7 temperature  values  for  each  entry.  This  allows  the  storage  of 
1620  items  of  data  specified  by  constituent  or  wavenumber,  with  seven 
temperature  dependent  values  for  each  item.  The  coefficient  data  may  be 
either  absorption  coefficient  or  line  density  (1/d)  data,  and  the  position 
for  the  storage  of  data  for  each  constituent  is  flexible.  The  position  of 
the  initial  coefficient  in  the  table  for  each  constituent  is  specified  to 
be  IC0EF  (I)  for  absorption  coefficients  of  IFIN  (I)  for  line  density.  The 
number  of  temperatures  at  which  data  is  stored  is  specified  as  NACT  (I)  for 
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absorption  coefficients  or  NFINT  (I)  for  fine  structure  parameters.  The 
coefficient  data  is  stored  in  the  block  data  statements,  so  the  problems 
associated  with  the  coefficient  storage  limitations  need  only  arise  if 
this  data  is  modified. 

The  maximum  number  of  spectral  intervals  which  may  be  used  is  440 
( (ENUF-ENUI )/DENlI  +1  <440).  If  this  limit  is  exceeded,  the  program  will 
stop  and  print  a diagnostic  message. 

The  maximum  number  of  combined  spectral  intervals  for  all  gaseous 
constituents  must  not  be  greater  than  840  (NUDIM).  If  this  limit  is  ex- 
ceeded in  M0DL3A,  the  calculation  will  stop  and  a diagnostic  message  will 
be  written.  The  limit  cannot  be  exceeded  in  M0DEL3  since  the  storage  is 
allocated  to  contain  the  calculation  within  the  desired  dimension.  How- 
ever, an  additional  limitation  is  placed  on  M0DEL3  that  the  number  of  line 
groups  (IRAD)  must  not  exceed  10.  If  it  does,  a diagnostic  message  is 
written  and  the  program  returns  to  try  to  read  the  next  case. 

Since  a line  of  sight  may  be  traversed  several  times  for  M0DEL3  in 
order  to  complete  the  spectral  range  desired  within  the  NUDIM  storage 
limitation,  the  spectral  output  will  only  be  for  the  spectral  range  of  the 
last  traverse  if  NSPEK  = 2 or  4.  In  addition,  the  line  of  sight  output 
using  NLINE  = 2 or  3 will  be  limited  to  first  traverse. 

The  restrictions  on  modifying  the  absorption  coefficient  data  are 
described  in  Item  15  of  SPECIAL  OPTIONS. 
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RESTART  PROCEDURE 

The  restart  procedure  will  depend  upon  the  type  of  output  desired. 

For  example,  if  spectral  output  for  the  entire  pnblem  is  desired,  the 
radiation  calculation  must  be  completed  at  one  time.  However,  in  the 
case  for  which  the  total  irradiance  is  of  primary  interest,  the  program 
may  be  restarted.  In  any  case,  a few  minutes  may  be  saved  if  tape  10  is 
saved  so  that  the  line  of  sight  property  data  does  not  have  to  be  recom- 
puted. Therefore,  instructions  should  be  given  to  save  tape  10  if  a 
lengthy  problem  is  being  run.  If  the  problem  is  to  be  restarted  from  a 
previously  generated  tape  10,  IFLOW  should  be  set  at  1,  and  only  the 
first  6 data  cards  are  required. 

The  irradiance  is  continually  summed  throughout  the  calculation  with 
output  depending  upon  the  output  option  (NLINE)  selected.  If  NLINE  = 0, 
the  accumulated  irradiance  (FLUX)  is  output  at  the  end  of  each  6 increment, 
but  with  all  the  other  NLINE  options  the  accumulated  irradiance  is  output 
after  each  line  of  sight.  When  the  program  is  restarted,  the  upper  limit 
of  the  final  value  of  THETA  which  is  completely  finished,  should  be  used 
as  the  initial  value  of  THETA  (THETAI)  when  the  problem  is  restarted. 

When  the  restarted  run  is  complete,  the  irradiance  (FLUX)  can  be  added 
to  that  predicted  by  the  first  run  to  obtain  the  total  irradiance  for  the 
problem. 
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DEFINITION  OF  TERMS 

MENMONIC 

MATH 

SYMBOL 

DEFINITION  UNITS 

ABSCL(I) 

Absorption  coefficient  scale  factor 
for  constituent  I. 

AC 

k(STP) 

Absorption  coefficient  at  standard  (l/cm)STP 

temperature  and  pressure 

ACMIN 

Minimum  absorption  coefficient  to  (l/cm)STP 

be  considered  in  radiation  calcula- 
tions. Input  on  card  6. 

ACO(I) 

Summation  of  the  collision  broadened 
fine  structure  parameter  (Yc/d)  times  ku. 

ADO(I) 

Summation  of  the  Doppler  broadened 
fine  structure  parameter  (Yg/d)  times  ku. 

AVGG 

Average  of  spectral  transmissivities 

\ 

I t(s  »w)/  ((Or  - 10-*)  J 

\Wi  / 

BCCO 

Vc 

Collision  half-width  for  carbon-  1/cm 

monoxide. 

BCC02 

\ 

Collision  half-width  for  carbon-  1/cm 

dioxide. 

BCDI 

ac 

Collision  broadened  fine  structure 
term  Yc/d. 

BCDICO 

ac 

Collision  broadened  fine  structure 
term,  y /d,  for  carbon  monoxide. 

BCH20 

Yc 

Collision  half-width  for  water  1/cm 

vapor. 

BD 

yd 

Doppler  half-width.  1/cm 

BDCO 

Doppler  half-width  for  carbon- 
monoxide  divided  by  to. 

BDC02 

Doppler  half-width  for  carbon- 
dioxide  divided  by  co. 

BDH20 

Doppler  half-width  for  water 
vapor  divided  by  to. 
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UNITS 


BDDI 

Doppler  broadened  fine  structure 
term  yD/d. 

BETA 

Flow  field  propellant  or  gas 
identification. 

C0EF(L,K) 

Absorption  coefficient  (k)  c'  line 
density  (1/d)  in  the  table  of 
coefficients  loaded  at  the  tempera- 
ture index  L for  table  position 
index  K. 

COLT 

Collision  broadening  term  in  the 
combined  optical  depth,  Y. 

CRATIO 

Radial  location  of  carbon  mole  fraction 
as  a fraction  of  the  miximum  radius. 
Input  on  card  14. 

CTCD 

Spacing  between  flow  field  axis  and 
the  Z-axis.  Input  on  card  19. 

cm  or  in. 

DASH 

Incremental  value  of  Planck's 
function  times  the  differential 
of  the  transmissivity. 

watts/cm-s  r 

DATE 

Date  of  run. 

DENU 

Aw 

Increment  for  the  spectral  summation 

1/cm 

DFO 

Incremental  value  of  the  optical 
depth,  D*. 

DIN 

1/d 

Line  density 

DINH20 

Fraction  of  1/d  used  as  l/d0  for 
water  vapor. 

DINSCL(I) 

Fine  structure  parameter  scale 
factor  for  consittuent  I. 

DINV 

1/d 

Line  density. 

cm 

DINVCO 

1/d 

Line  density  for  carbon 
monoxide 

cm 

DOPT 

Doppler  term  in  the  combined 

optical  depth,  Y. 
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DEFINITION 

UNITS 

DM0 

X 

Sum  of  the  optical  depths  for  all 
constituents. 

DPHI 

A<p 

Integration  increment  for  the  4> 
integral 

deg 

DPHIS 

Storage  for  input  value  of  DPHI. 

deg 

DS 

Integration  interval  along  the 
line  of  sight 

in  or  cm 

DSFACT 

Increment  of  SFACT  and  TSFACT  for 
a single  line  of  sight. 

(sin  6 cos  0 A0A4>/n) 

DSFLUX 

Incremental  value  of  the  radiation 
flux  for  a geometric  increment 
along  a line  of  sight 
(NJAt  sin  0 cos  0A0A4>Aoj) . 

watts/cnf 

DSS 

As 

Geometric  integration  step  alonq  the 
line  of  sight.  This  differs  from  DS 
in  that  it  may  be  a sum  of  individual 
geometric  values  over  a temperature 
increment. 

DTHETA 

AO 

Increment  for  the  summation  over  0 

deg 

DUCARB 

Incremental  value  summed  to  give 
UCARB. 

gm/cm2 

DUCO 

Incremental  value  summed  to  give 
UCO. 

atm-cm 

DUC02 

Incremental  value  summed  to  give 
UC02. 

atm-cm 

DUH20 

Incremental  values  summed  to  aive 
UH20. 

atm-cm 

ELAM 

Wavelength 

mi  cron 

END 

b) 

Wavenumber  at  the  center  of  a Aw 
wavenumber  interval. 

1/cm 

ENUF 

Upper  limit  of  the  w integral. 

1/cm 

enui 

"Wi 

Lower  limit  of  the  w integral. 

1/cm 
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EPS 

e 

Parameter  used  for  gn  in 

computing  l/dn  for  carbon 
di oxi de- 

ETA(I.K) 

value  of  n in  the  K-position 
on  the  Ith  Z-plane. 

ETAC(I.K) 

Angular  position  of  the  Kth 
carbon  value  in  the  Ith 
Z-plane.  Input  on  card  13. 

EXP0(1,N) 

Fraction  of  the  water  vapor 
absorption  coefficient  for  the 
nth  line  group. 

EX?0(2,N) 

Fraction  of  the  carbon  dioxide 
absorption  coefficient  for  the 
nth  line  group. 

EXP0(3,N) 

Fraction  of  the  carbon  monoxide 
absorption  coefficient  for  the 
nth  line  group. 

EYE 

N°„ 

Plack's  function 

FISCAL 

Scale  Factor  for  the  first  mole 
fraction. 

F2SCAL 

Scale  factor  for  the  second  mole 
fraction. 

F3SCAL 

Scale  factor  for  the  third  mole 
fraction. 

F4SCAL 

Scale  factor  for  the  fourth  mole 
fraction. 

F5SCAL 

Scale  factor  for  the  fifth  mole 
fraction. 

F6SCAL 

Scale  factor  for  the  sixth  mole 
fraction. 

FC(I,L,K) 

Carbon  mole  fraction.  Input  on 
Card  14. 

FF(J.L) 

Mole  fraction  of  the  Lth  constituent 
being  read  in  from  the  flow  field 
tape  at  radial  point  J. 

RTR  013-1 
UNITS 


deg 

deg 


watts/cm2 -sr 


L 
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MATH 

SYMBOL 

DEFINITION 

UNITS 

FLUFF 

Product  of  integration  intervals 

(A0A$A to)  . 

FLUX 

q/A 

Radiant  heat  transfer  rate  per 
unit  area.  When  used  as 
intermediate  output.it  is  the 
summation  of  the  radiant  heat 
transfer  rate  up  to  the  point 
at  which  it  is  printed. 

watts/cm 

F0(I) 

X* 

Optical  depth  for  the  weak 
line  approximation  at  the 
Ith  wavenumber  interval. 

FOC(I) 

X* 

Optical  depth  for  carbon  particles 
at  the  Ith  wavenumber  interval. 

FW(L) 

Mole  fraction  of  the  Lth  constituent  at 
line-of-sight  point  under  consideration. 
(Index  based  on  input  order  on  card  9). 

the 

FWSUM(II) 

s 

Sum  of  mole  fractions  for  the  I Ith 
constituent  over  a temperature 
Increment. 

G1 

thru 

G6 

^I.j^STP 

Coll  is ion  broadening  term  for 
the  unidentified  gas  represented 
by  PREM.  Contained  in  program  data 
statement.  Can  be  modified  using 
card  25. 

1/cm-atm 

GAM 

Ratio  of  temperature  to  standard 
temperature  (273K)  to  the  1/2  power. 

GNEW(I) 

t(s  ,6l>) 

Transmissivity  for  the  position  s and 
the  Ith  wavenumber  interval  centered 
at  co. 

HANG 

Flow  field  sector  angle.  Input  on 
card  19. 

HDG 

Flow  field  title. 

I BLOCK 

Number  of  blocking  circles.  Input  on 
card  2. 

ICARB 

Carbon  option.  Input  on  card  2. 

ICOEF(L) 

Position  of  the  first  absorption 
coefficient  for  the  constituent  L 
in  the  table  of  coefficients. 
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ICOEFF 

Option  to  modify  absorption 
coefficient  data.  Input  on  card  2. 

IDENTC(I) 

Constituent  names  recognized  by  the 
program.  The  order  of  IDENTC  list 
in  the  MAIND  deck  governs  the  order 
of  operation  of  the  program. 

IDENU 

Aw 

Increment  for  the  summation  over  w 

1/cm 

IDNU(I) 

Wavenumber  interval  in  the 

table  of  coefficients  for  the  Ith 

constituent. 

1/cm 

IENUI 

wi 

Lower  limit  of  the  w integration 

1/cm 

IENUF 

“f 

Upper  limit  of  the  w integration. 

1/cm 

IFLOTP 

Index  to  control  production 
of  a flow  field  tape.  Input 
on  card  2. 

I FLOW 

Index  to  specify  manner  of 
flow  field  input.  Input  on 
card  2. 

IMOD 

Number  of  constituents  for 
which  coefficient  data  is  to  be 
modified.  Input  on  card  26. 

I PLANE 

Option  for  specification  of 
plane  orientation.  Input  on  card 
2. 

I RAD 

Number  of  line  groups  to  be 
used  in  Model  3 radiation 
calculation.  Input  on  card  2. 

ISCALA 

Absorption  coefficient  scaling 
option.  Input  on  card  2. 

ISCALF 

Geometry  and  gas  property  scaling 
option.  Input  on  card  2. 

ISHOCK 

Shock  plane  option.  Input  on 
card  2. 

ISYMSC 

Three  dimensional  flow  field  parameter. 
Input  on  card  19. 

INUITS 

Input  units  option.  Input  on  card  2. 
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Hij 


I 


I 


IZF 

IZI 

IZSH 


IZSUM 

JPLOT(I) 

JFK 

KASE 

KCARB 

KIND  ( I ) 
KOM(I) 


LBLOCK 


LIBCON 


Index  of  final  Z-plane  to  be 
considered  in  an  input  flow  field. 

Input  on  card  2. 

Index  of  initial  Z-plane  to  be 
considered  in  the  input  flow  field. 

Input  on  card  2. 

Index  of  the  initial  Z-plane  in 
the  flow  field  for  which  shock 
data  is  specified.  Input  on 
card  15. 

Index  of  the  most  recently  loaded 
(input)  Z plane  referenced  to  the  complete 
gas  flow  field. 

Plotting  control  input  on  card  32. 

Index  used  to  indicate  the  number  of 
geometric  steps  within  a temperature 
step  interval. 

Case  number.  Input  on  card  2. 

Index  of  the  carbon  mole  fraction 
in  the  flow  field  storage  when 
ICARB>0. 

Indicates  if  the  Ith  blocking  circle 
is  a disk  (0)  or  a hole  (1 ). 

Constituent  name.  Each  name  consists 
of  two  A6  words,  but  none  of  the  names 
use  more  than  7 characters.  The  name 
cannot  be  changed  unless  the  program 
used  to  generate  the  flow  field  input 
tapes  is  modified.  Input  on  card  9. 

Index  to  indicate  if  the  line  of 
sight  has  been  blocked  by  a blocking 
circle.  If  it  has  been  blocked, 

LBLOCK  will  be  the  index  of  the  block- 
ing circle  which  blocked  it. 

Number  of  constituents  identified 
in  the  program. 
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LOSPT 

LOSFLX 

LPT 

LSET 

NABS 

NACT(I) 

NET 

NETA(I) 

NETAS 

NFINT(I) 

NFLOW 

NLINE 

NOCON 

NOETAS(I) 


f-IATH 

SYMBOL  DEFINITION 


Index  for  the  position  of  points 
along  a line  of  sight. 

Summation  of  the  heat  flux  per  unit 
area  along  a line  of  sight  (EN^At  sine 
cos  6 A6A(}>  Aco) 

Number  of  line  of  sight  points  on  a 
tape  record. 

This  index  provides  the  number  order 
for  the  lines  of  sight  on  the  line  of 
sight  tape.  A zero  value  indicates  the 
end  of  all  line  of  sight  records  on  the 
tape,  and  a negative  value  indicates  that 
the  line  of  sight  has  either  been  blocked 
or  missed  the  flow  field. 

Absorption  coefficient  output  option. 
Input  on  card  2. 

Number  of  temperature  values  at  which 
absorption  coefficients  are  stored  for 
the  Ith  constituent. 

Number  of  n values  in  a Z-plane. 

Number  of  n values  in  the  Ith  Z-plane. 

Number  of  n values  in  a Z-plane. 

Number  of  temperature  values  at 
which  fine  structure  parameters 
are  stored  for  the  Ith  constituent. 

Flow  field  output  option. 

Input  on  card  2. 

Radiation  output  option  index.  Input 
on  card  2. 

Number  of  constituents  stored  on  flow 
field  Input  tape. 

Number  of  n values  input  at  the  Ith 
Z-plane  for  shock  data.  Input  on 
card  15. 


UNITS 


watts/cm2 
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MATH 

SYMBOL 

DEFINITION 

UNITS 

NOETAC(I) 

Number  of  n values  used  for  loading  carbon 
mole  fraction  at  the  Ith  carbon  Z-plane. 
Input  on  card  12. 

NOKON 

Number  of  constituents  to  be  considered 
for  the  radiation  calculation.  Input 
on  card  2. 

NOPS 

Number  of  points  on  a radial  line  in 
the  flow  field. 

N0PSFC(I,L) 

Number  of  points  at  which  carbon 
mole  fraction  is  specified  at  the  Ith 
carbon  Z-plane  and  Lth  n position. 
Input  on  card  13. 

NPLOT 

Number  of  plots  desired.  Input  on 
card  2. 

NR(I) 

Number  of  radial  points  at  ZZ ( I ) . 

NR(I ,J) 

Number  of  radial  points  at  ETA(I,J). 

NSPEK 

Spectral  output  option.  Input 
on  card  2. 

NU 

to 

Wavenumber  at  the  center  of 
a Acj  wavenumber  interval. 

1/cm 

NUDIM 

Number  of  storage  positions  reserved 
for  optical  depth  parameters  of 
gaseous  radiators. 

NUF 

w 

Upper  limit  of  the  w integral. 

1/cm 

NUI 

to 

Lower  limit  of  the  u integral. 

1/cm 

NUL(I) 

Value  of  the  minimum  wavenumber 
in  the  table  of  the  absorption 
coefficients  for  constituent  I. 

NUMBER(I) 

Position  of  the  constituent  in 

the  IDENTC  list  for  which  the 
coefficient  data  is  to  be 
modified,  input  on  card  26. 


49 


REMTECH  INCORPORATED 


RTR  013-1 


MATH 

MNEMONIC  SYMBOL  DEFINITION 


UNITS 


NUMCON(I) 


NUU(I) 

NFZS 

NZFC 

NZSH 

0D(I) 

X 

OMEGA 

P(I.L) 

P 

P(IXJ) 

P 

PC 

PCO 

Pi 

PC02 

Pi 

PCORR 

PH2 

Pi 

Index  of  the  position  of  the  Ith 
constituents  in  the  identifying 
data  list  (IDENTC)  in  the  program. 

Value  of  the  maximum  wavenumber  in 
the  absorption  coefficient  table  for 
constituent  I. 


Index  of  the  last  Z-plane  in  the  flow 
field  for  which  shock  data  is  specified. 

Number  of  Z-planes  used  for  adding  carbon 
mole  fraction.  Input  on  card  11. 

Number  of  Z-planes  for  which  shock  data 
is  loaded.  Input  on  card  15. 

Optical  depth  for  the  index  I.  The 
I index  used  in  the  spectral  summation 
is  a function  of  both  wavenumber  and 
constituent. 


Inclination  angle  of  the  plane  of 
interest. 

Pressure  at  the  Lth  point  on  ZZ(I). 

Pressure  at  the  Jth  point  on 
ETA( I , K) . 

Carbon  partial  pressure  corrected 
for  temperature.  (Note  that  carbon 
particles  are  trested  as  a gas  of 
molecular  weight  12.011) 

Carbon  monoxide  partial  pressure 
corrected  for  temperature. 

Carbon  dioxide  partial  pressure 
corrected  for  temperature. 

Pressure  corrected  to  standard 
temperature.  This  is  actually 
a ratio  of  the  gas  density  to  the 
density  at  standard  conditions. 

Hydrogen  partial  pressure  corrected 
for  temperature. 


deg 

psfa  or  atm 
psfa  or  atm 

atm 

atm 

atm 

atm 

atm 


i 
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UNITS 

PH20 

Pi 

Water  vapor  partial  pressure 
corrected  for  temperature 

atm 

PHI 

Angle  between  the  positive  V axis  and 
projection  of  a line  of  sight  in  the 
UV  plane. 

deg 

PHIF 

Upper  limit  of  the  <j>  summation. 

deg 

PHIFS 

Storage  for  input  value  of 
PHIF. 

deg 

PHI  I 

Lower  limit  of  the  <J>  summation. 

deg 

PHI  I S 

Storage  for  input  value  of  PHI I. 

deg 

PN2 

pi 

Nitrogen  partial  pressure 
corrected  for  temperature. 

atm 

P02 

pi 

Oxygen  partial  pressure 
corrected  for  temperature. 

atm 

PREM 

Partial  pressure  of  all  remaining 
gases  not  accounted  for  by  the 
identified  constituents. 

atm 

PSCAL 

Pressure  scale  factor. 

PSFINV 

Inverse  of  the  integral  of 
sine  cose  de  d<p. 

PW 

P 

Pressure  at  the  line  of  sight 
point  under  consideration. 

atm 

PWSUM 

Sum  of  the  pressures  over  a 
temperature  increment. 

atm 

R(I.L) 

Radius  of  the  Lth  point  on  ZZ ( I ) . 

in/cm 

R(I»K,J) 

Radius  of  the  Oth  point  on  ETA(I,K). 

in/cm 

RADAVG 

Averaae  radiance  over  all  lines  of 
sightT  (FLUX/(TSFACT*ir)) 

watts/cm2-sr 

RADLOS 

Radiance  for  a particular  line 
of  sight.  (SUM  * Aw) 

watts/cm2-sr 

I 
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UNITS 

RADPLM 

Average  radiance  of  the  lines  of 
sight  which  went  through  the  gas 
flow  field.  (FLUX/(SFACT*tt)) 

watts/cm2-sr 

RAM 

Ratio  of  standard  temperature  (273K) 
to  the  gas  temperature  to  the 
1/2  power. 

RC(I) 

Radius  of  the  Ith  blocking  circle. 

in  or  cm 

RHOC 

P 

Carbon  particle  density. 

gm/cm3 

RSCAL 

Geometric  scale  factor  for  problem 
dimensions. 

RSHOCK  (1,0) 

Shock  radial  position.  Input 
on  card  16. 

cm  or  in 

S 

s 

Distance  along  a line  of  sight  from 
the  point  of  interest. 

cm 

SANKA 

Product  of  the  angle  factors  and 
integration  intervals  (sine  cos0 
A0A$Aw) . 

SFACT 

Summation  of  cose  sin©  A0  A^/ir 
for  the  lines  of  sight  which  intersect 
the  exhaust  plume. 

SFLXTH(I) 

Spectral  flux  in  the  Ith  wave 
number  increment  for  the  current 
value  of  0. 

watts/cm2-cm 

SFLXTO(I) 

Sum  of  the  spectral  flux  in  the  Ith 
wavenumber  increment.  Also  used  as 
storage  for  spectral  output. 

watts/cm2-cm' 

SINKO 

The  product  of  cose  sine. 

SMAX 

L 

Input  value  of  the  upper  limit  of 
the  s summation.  Input  on  card  6. 

in  or  cm 

SMIN 

L 

Value  of  the  upper  limit  of  the  s 
summation  after  the  line  of  sight 
has  been  checked  for  obstruction. 

cm 
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SUM 

Summation  over  a line  of  sight  of  the 
product  of  the  Planck  function  (N°w) 
and  the  differential  transmissivity 
(At). 

watts/cm2 -s 

SUMOLD 

Storage  location  for  the  radiation 
flux  thru  the  previous  line  of  sight 
step. 

SZERO 

Lower  limit  of  the  summation  over  s. 
In  this  program,  it  is  set  at  zero. 

T(I,L) 

T 

Temperature  of  the  Lth  point  on  ZZ ( I ) . 

R or  K 

T(I.M) 

T 

Temperature  of  the  Jth  point  on 
ETA(I,K). 

R or  K 

TANGLE 

Input  value  of  (sinecose  A0A<f>)  to 
be  used  in  computing  DPHI. 

sr 

TDIFF 

Temperature  change  limit  used  in 
determining  integration  interval 
along  a line  of  sight.  Input  on 
card  6. 

K or  R 

TEMAC ( I , K) 

K values  of  temperatures  at  which 
absorption  coefficients  for  the 
constituent  I are  listed  in  the 
coefficient  table. 

K 

TEMFIN(I.J) 

0 values  of  temperature  at  which 
line  density  (1/d)  for  a constituent 

1 are  listed  in  the  coefficient  table. 

K 

THETA 

6 

Angle  between  the  surface  normal  at 
the  point  of  interest  and  a line  of 
sight  ray  from  the  point  of  interest. 

deg 

THETAI 

0i 

Lower  limit  of  the  0 summation 

deg 

THETAF 

6f 

Upper  limit  of  the  0 summation 

deg 
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TINDEX 

Index  used  for  temperature  step 
computation  to  determine  the  type 
of  steps  being  considered:  minus  one 
indicates  the  initial  step  on  the  line 
of  sight;  zero  indicates  the  initial 
step  in  a new  temperature  increment;  and 
one  indicates  a step  within  a temperature 
increment  that  is  being  summed. 

TSCAL 

Temperature  scale  factor. 

TSFACT 

Summation  of  cos 9 sine  A9  A<J>/tt 
for  all  lines  of  sight. 

TRUESF 

Exact  value  for  the  integral  for  sine 
cos0  d9  d<J>/ir  for  all  lines  of  sight. 

TWSTO 

Temperature  at  the  initial  geometric 
step  on  a temperature  step. 

K 

TWSUM 

Sum  of  temperatures  over  a temperature 
increment. 

K 

TW 

T 

Temperature  at  the  line  of  sight  point 
under  consideration. 

K 

UCARB 

u 

Density-path  length  product  for 
carbon  particles. 

gm/sq-cm 

UCO 

u 

Pressure-path  length  product  for 
carbon  monoxide. 

(cm-STP) 

UC02 

u 

Pressure-path  length  product  for 
carbon  dioxide. 

(cm-STP) 

UH20 

u 

Pressure-path  length  product  for 
water  vapor. 

(cm-STP) 

ULX 

Direction  cosine  of  U-axis  with 
respect  to  X-axis. 

ULY 

Direction  cosine  of  U-axis  with 
respect  to  Y-axis. 

ULZ 

Direction  cosine  of  U-axis  with 
respect  to  Z-axis. 

VLX 

Direction  cosine  of  V-axis  with 
resp  ct  to  X-axis. 

VLY 

Direction  cosine  of  V-axis  with 
respect  to  Y-axis. 
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VLZ 

Direction  cosine  of  the  V-axis  with 
respect  to  the  Z-axis. 

WLX 

Direction  cosine  of  the  W-axis  with 
respect  to  the  X-axis. 

WLY 

Direction  cosine  of  the  W-axis  with 
respect  to  the  Z-axis. 

WLZ 

Direction  cosine  of  the  W-axis  with 
respect  to  the  Z-axis. 

X 

X 

X coordinate  at  the  line  of  sight 
point  under  consideration. 

cm 

xc(i) 

X coordinate  of  the  Ith  blocking 
circle. 

in  or  cm 

XP 

X coordinate  of  the  point  of  interest. 

in  or  cm 

Y 

Y 

Y coordinate  at  the  line  of  sight 
point  under  consideration. 

cm 

YC  ( I ) 

Y coordinate  of  the  Ith  blocking 
circle. 

in  or  cm 

YP 

Y coordinate  of  the  point  of  interest. 

in  or  cm 

Z 

Z 

Z coordinate  at  the  line  of  sight 

cm 

ZC(I) 

Z coordinate  of  the  Ith  blocking 
circle. 

in  or  cm 

ZP 

Z coordinate  of  the  point  of  interest. 

in  or  cm 

ZZ(I) 

Z coordinate  of  the  Ith  Z plane. 

in  or  cm 

ZZFC(I) 

Z coordinate  of  the  Ith 

carbon  Z plane.  Input  on  card  12. 

in  or  cm 
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Card 

Format 

Col 

Computer 

Value 

Description 

Name 

1 

8A10 

1-80 

TITLE 

Problem  identification 

2 

16,2013 

1-6 

KASE 

<0 

Flowfield  input  from  cards 

0 

See  Note  B on  page  59  . 

>0 

Case  number  of  flow  field  on  tape  8 
If  IFLOTP  >0  use  tape  11 

8 

IFLOTP 

0 

Flowfield  not  modified 

Flowfield  to  be  modified  as  required 
and  output  on  tape  if  IFL0TP>0 

1 

First  case  on  tape  8 

2 

Other  than  first  case  on  tape  8 

9 

I FLOW 

<1 

Use  previous  LOS  tape  mounted  on 
tape  10.  See  notes  on  following  page 

2 

Axisymmetric  flow  field  input. 

3 

3-Dimensional  flow  field  input. 

12 

N0K0N 

Number  of  constituents  to  be  stored 
in  flow  field  (maximum  of  6). 

15 

IPLANE 

0 

Read  OMEGA  for  plane  of  inclination. 

1 

Read  direction  cosines  of  U,  V, 
W,  axes. 

17,18 

IBLOCK 

Number  of  blocking  circles  used 
(50  maximum). 

21 

ICARB 

0 

Read  carbon  mole  fraction  from  the 

flow  field  input  (if  carbon  is  called 
for  and  available  in  flow  field). 

1 Use  constant  carbon  mole  fraction 
FC=constant 

2 Use  carbon  mole  fraction  as  a 
function  of  radius:  FC=f(R)  for 
axi symmetric  flowfields  and  FC=f 
(n,R)  for  3-D  flowfields. 

3 Use  carbon  mole  fraction  as  a 
function  of  Z and  R:  FC=f(Z,R)  for 
axi symmetric  flow  fields  and  FC=f 
(Z,n,R)  for  3-D  flow  fields. 
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Card  Format  Col  Computer  Value 
Name  

2 

(cont.)  24  NLINE  0 

1 

2 

3 

26  NSPEK  0 

1 

27  0 

1 

2 

3 

4 

29,30  NPLOT 
33  NFLOW  0 

1 

36  NABS  0 

1 

39  ISHOCK  -1 

0 

1 

42  I UNITS  0 

1 


Description 


Output  cumulative  flux  and  shape 
factor  at  each  6 increment. 

Output  cumulative  flux,  shape 
factor,  LOS  flux,  and  LOS  intensity 
for  each  LOS. 

Same  as  1 plus  properties  at  each 
radiation  step. 

Same  as  2 plus  properties  at  each 
geometric  step. 

Spectral  output  per  unit  wavelength 
Spectral  output  per  unit  wavenumber 

No  spectral  output  desired 
Spectral  radiance  per  steradian  - 
short  form 

Spectral  radiance  per  steradian  - 
long  form 

Spectral  radiarce  - short  form 
Spectral  radiance  - long  form 

Number  of  plots  desired. 

Short  form  of  flow  field  output  with 
only  first  and  last  radius  points. 
Complete  flow  field  output. 

Short  form  of  absorption  coefficient  out- 
put with  only  first  and  last  values. 
Complete  absorption  coefficient  output. 

Shock  location  input  with  the  flowfield. 
Shock  plane  not  specified. 

Shock  plane  data  input. 

English  units  (inches,  °R,  lbf/ft2) 
used  in  all  input  and  output  data 
except  absorption  coefficients. 

MKS  units  used  in  all  input  and  output 
data  (cm,  °K,  atm). 
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Card  Format 

Col 

Computer 

Value 

Description 

Name 

2 

45 

ISCALF 

0 

No  flow  field  scale  factors  used. 

(cont. ) 

1 

Flow  field  scale  factors  input. 

47,48 

IZI 

0 

Read  first  Z-plane  in  the  flow  field 

> 0 

Index  of  first  Z-plane  to  be  read. 

50,51 

IZF 

0 

Read  to  end  of  flow  field. 

> 0 

Index  of  last  Z-plane  to  be  read. 

53,54 

I RAD 

0 

Use  Model  3A  radiation  calculation 

> 1 

Use  Model  3 radiation  calculation 
with  line  groups  n = IRAD 

.. 

(maximum  of  10) 

57 

IC0EFF 

0 

Use  absorption  coefficient  data 
in  program. 

1 

Input  modified  absorption 
coefficient  data. 

-1 

Use  analytical  expression  for  H,0 
1/d. 

60 

ISCALA 

0 

No  absorption  coefficient  scale 
factors  used 

1 

Absorption  coefficient  scale  factors 
input. 

NOTE: 

Procedure  for 

using  Previously 

Generated 

Tape  10 

A.  If  a change  is  desired  in  problem  limits: 

1.  KASE  must  match  the  absolute  value  of  the  case  number  used  in 
generating  tape  10. 

2.  Cards  3-6  are  used  to  specify  new  limits.  (Increments  must  be  the  same) 
The  values  used  must  produce  line-of-sight  angles  which  exist 

on  the  tape. 

3.  Omit  cards  7-16  and  18-22.  Remaining  cards  are  used  as 
required  by  options  chosen. 

B.  If  no  change  is  desired  in  problem  limits: 

1.  Set  KASE  = 0. 

2.  Omit  cards  3-16  and  18-22.  Remaining  cards  are  used  as  required 
by  options  chosen. 
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Card 

Format 

Col 

Computer 

Name 

Units 

3 

3E10.6 

1-10 

THETA I 

Deg 

11-20 

THETAF 

Deg 

21-30 

DTHETA 

Deg 

4 

3E10.6 

1-10 

PHII 

Deg 

11-20 

PHIF 

Deg 

21-30 

DPHI 

Deg 

31-40 

TANGLE 

Sr 

5 

3E10.6  1-10 

ENUI 

_ l 

cm  x 

11-20 

ENUF 

cm" 

21-30 

DENU 

cm"1 

Description 


Initial  value  = 8^ 

Final  Value  = 6f 
Increment  = A8  T 

_ . . , . j Note:  <j>  is  clockwise 

Initial  value  = ^ from  u when  viewed  in  the 
Final  value  = positive  W direction. 
Increment  = A <P 

If  this  value  is  specified,  A<p  will 
be  adjusted  so  that  cose  sin8A0A<J)  = 

TANGLE. 

Wavenumber  lower  limit, 

Wavenumber  upper  limit,  wf 
Increment,  Aw.  Must  be  an  interger  multiple 
of  the  increment  in  the  absorption 
coefficient  table. 


6 


3E10.6  1-10  TDIFF 

11-20  SMAX 


R or  K Temperature  step  size  for  integration  on  S. 
in  or  cm  Upper  limit  on  LOS  distance,  L. 


21-30  ACMIN  (cm"1)  Minimum  absorption  coefficient  to  be  used 

b p in  radiation  calculation. 


7 6E10.6  1-10  XP  in  or  cm  X) 

11-20  YP  in  or  cm  Y/  Coordinates  of  the  point  of  interest. 

21-30  ZP  in  or  cm  Z* 

31-40  DS  in  or  cm  Geometric  increment  along  LOS,  As. 

NOTE:  If  I PLANE  = 0,  read  8,  but  if  IPLANE  = 1,  read  8A 


8 

E10.6 

1-10 

OMEGA  Deg 

Inclination  of  the  plane  of  interest,  ft 

8A 

9E8.4 

1-8 

ULX 

Direction  cosines 

for  U-X  axes. 

9-16 

ULY 

i 

i 

U-Y 

17-24 

ULZ 

U-Z 

25-32 

VLX 

V-X 

33-40 

VLY 

V-Y 

41-48 

VLZ 

V-Z 

29-56 

WLX 

W-X 

57-64 

WLY 

r 

W-Y 

65-72 

WLZ 

Di rection 

cosines 

for  W-Z  axes. 

9 

12A6 

1-12 

K0N  (1 ,2) 

First  constituent 

Note:  The  constituent 

13-24 

K0N (3,4) 

Second  constituent 

names  must  be  left  justi- 

25-36 

K0N  5,6) 

Third  constituent 

fied  and  must  be  in  the 

37-48 

K0N(7,6) 

Fourth  constituent 

following  forms:  H201 (G ) , 

49-60 

K0N(9,10) 

Fifth  constituent 

Cl  02(G), C101 (G),C1(S), 

61-72 q 

KON  (11 ,12) 

Sixth  constituent 

N2(G),  H2(G),  02(G), 

HCL,  HF. 
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INPUT  CARDS  (Continued) 


Card  Format  Col  Computer  Units 
Name  

USE  CARD(S)  10  IF  IBL0CK>0.  1 = 1,  IBLOCK 

10  4E10.6 , 1-10  XC(I)  in  or  cm 


4E10.6 , 1-10  XC(I) 

15  11-20  YC(I) 

21-30  ZC(I) 
31-40  RC( I ) 
45  KIND ( I ) 


Description 


in  or  cm  Xj 

in  or  cm  Y;  Coordinates  of  the  center  of  the 
in  or  cm  Z;  blocking  circle 
in  or  cm  Blocking  circle  radius 

Type  of  blockage:  0 = disk;  1 = hole 


USE  CARDS  11  THRU  14  FOR  ADDING  CARBON  MOLE  FRACTION  IF  ICARB  >0. 
The  order  of  using  these  cards  is  as  follows: 


I FLOW 

ICARB 

Procedure 

2 

: or  3 

1 

Read  card  14  one 

time  for  FC  (1 , 1 , 1 ) 

2 

2 

Read  card  13  one  time  for  NOPSFC  (1,1) 
Read  card(s)  14, K = 1,  NOPSFC  (1,1) 

3 

Read  card  11  for  NZFC 

Read  cards  12  thru  14  for  I * 1,  NZFC  and  L = 1 

3 

2 

Read  card  12  for  NOETAC  (1) 

Read  card  13  and  14  for  I = 1 and  L = 1,  NOETAC  (1) 

3 

Read  card  11  for  NZFC 

Read  cards  12  thru  14  for  1=1,  NZFC 

Read  cards  13  and  14  for  1=1,  NZFC  and  L = 1,  NOETAC  (I) 

Card 

Format 

Col 

Computer  Units 

Name 

Description 

11 

12 

1,2 

NZFC 

Number  of  Z-planes  at  which  carbon  mole 
fraction  is  input  (maximum  of  10) 

12 

I2.E18.0 

1,2 

3-20 

NOETAC ( I ) 

ZZFC(I)  in  or  cm 

Number  of  n-planes  with  carbon  data  at 
this  Z-plane  (maximum  of  7).  Not 
required  for  IFLOW  = 2. 

Z coordinate  of  this  Z-plane 

13 

I2.E18.0 

1,2 

3-20 

NOPSFC ( I ,L ) 

ETAC ( I ,L ) Deg 

Number  of  radius  points  at  which  the  carbon 
mole  fraction  is  specified. 

Angle  of  this  ty-plane  (not  required 
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INPUT  CARDS  (Continued) 


1 


Card  Format  Col 

Computer  Units 

Name 

14  8E10.6  1-10 

FC(I ,L,1 ) 

11-20 

CRATI0 ( I ,L  ,1 ) 

21-30 

FC(I ,L,2) 

31-40 

CRATI0 ( I ,L ,2 ) 

41-50 

51-60 

61-70 

FC(I ,L,K) 

71-80 

CRATI0(I ,L,K) 

USE  CARDS  15  AMD  16  ONLY  IF  ISHOCK  = 1 


15 

2213 

1-3 

NZSH 

4-6 

IZSH 

7-9 

NOETAS(I) 

64-66 

NOETAS(NZSH) 

16 

8E10.6 

1-10 

11-20 

RSH0CK(1 ,1 ) in  or 

t 

71-80 

▼ 

RSH0CK(I,J) 

USE  CARD  17  ONLY  IF  ISCALF  = 1 


17  9E8.5  1-8  RSCAL 

9-16  TSCAL 
17-24  PSCAL 
25-32  FISCAL 
33-40  F2SCAL 
41-48  F3SCAL 
49-56  F4SCAL 
57-64  F5SCAL 
65-72  F6SCAL 


cm 


Description 


Carbon  mole  fraction  at  radius  point  1. 
Radius  ratio  (R/RMAX)  at  radius  point  1. 


i 


Repeat  with  as  many  cards  as  necessary 
up  to  a maximum  of  K=10. 


Number  of  consecutive  Z-planes  with 
shock  points  (maximum  of  20) 

Index  of  first  Z-plane  for  which  a shock 
is  specified 

For  I FLOW  = 3,  this  specifies  the  number 
of  n-planes  at  which  shock  points  are 
specified.  This  should  correspond  to 
the  n-planes  in  the  flow  field  at  the 
corresponding  Z-plane. 

Shock  radii-  list  for  reading  in  a double 
loop  going  thru  J = 1,  NOETAS(I)  first 
(inside)  and  1=1,  NZSH  second  (outside). 
For  IFLOW  = 2,  J = 1 only.  If  no  shock 
exists  at  a particular  value  of  n,  use 
RSHOCK  = 0. 


Dimension  scale  factor 

Temperature  scale  factor 

Pressure  scale  factor 

Mole  fraction  scale  for  constituent  1 

2 

3 

4 

5 

▼ 6 


USE  CARDS  18  THRU  22  ONLY  IF  KASE  <0  (FLOW  FIELD  INPUT  ON  CARDS) 


18  72H  1-72  HDG  Flow  field  description 

24H,I1,  1-24  BETA  Flow  field  gas  identification 

5X,2E10.6 


19 
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INPUT  CARDS  (Continued) 

Card  Format  Col  Computer  Units 

Name 


Description 


25  ISYMSC 


31-40  HANG 
41-50  CTCD 


0 Flow  field  of  center  engine  assumed 
to  be  symmetric  with  that  of  the  out- 
board engine  about  a plane  which  bisects, 
and  is  normal  to,  a plane  connecting 

the  engine  centerlines.  See  Special 
Option  21. 

1 No  special  assumptions  about  the 
center  engine  flow  field  or  no 
center  engine. 

Deg  Included  angle  of  flow  field  sector,  tfi. 
in  or  cm  Distance  from  flow  field  axis  to  Z-axis. 
For  a symmetrical  engine  arrangement, 
the  Z-axis  must  be  at  the  center  of 
symmetry. 


Repeat  cards  20  thru  22  as  required.  END  FLOW  FIELD  BY  A BLANK  CARP. 
20  E10.6  1-10  ZZ( I ) in  or  cm  Coordinate  of  this  Z-plarn 


20  El  0.6 

110,15 


18,20 


NR  ( I ) 
NETA(I) 


23,25  LS(I) 


E10.6 
HO, 15 


1-10  ETA(I.J) 


18-20  NR 
23-25  LS 


in  or  cm  Coordinate  of  this  Z-plane. 

I FLOW  = 2 only.  Number  of  radius  points 

I FLOW  = 3 only.  Number  of  q-planes 

Used  for  IFL0W  = 2 and  ISH0CK  =-l  only. 
Index  of  the  radial  point  just  inside 
of  a shock  position.  In  order  for  this 
option  to  allow  the  shock  location  to 
be  recognized,  LS(I)>2,  and  there  must  be 
two  points  on  the  outside  of  the  shock. 


USE  CARD  21  FOR  I FLOW  = 3 ONLY 

Repeat  cards  21  and  22  for  J = 1,  NOETA(I) 


Angle  of  the  J th  n-plane.  J=  1, 
NOETA(I) 

Number  of  radius  points 
ISH0CK  = -1  only.  Index  (L)  of  radial 
point  just  inside  the  shock  surface.  In 
order  this  option  to  perform  properly, 
there  must  be  at  least  two  radial  points 
on  each  side  of  the  shock. 


Repeat  card  22  L = 1,  NR(I,J) 
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INPUT  CARDS  (Continued) 


Card  Format 


Col 


Computer 

Name 


Units 


Description 


22 


9E8.5 


1-8 

9-16 

17-24 


I.J.L 

I.J.L 


in  or  cm  Radius  (J=l  for  IFL0W=2) 
R or  K Temperature 


P(l,J,L)  psfa  or  atm  Pressure 


25-32 

FI 

[I.J.L) 

First  \ 

33-40 

F2 

[I.J.L)  - 

Second 

41-48 

F3 

:i,j,l  - 

Third 

49-56 

F4 

[i.j.l) 

Fourth/ 

57-64 

F5 

I.J.L)  - 

Fifth  1 

65-72 

F6 

[I.J.L) 

Sixth  | 

Constituent  mole  fractions  - 
the  order  must  match  the 
order  used  on  card  9 for 
the  constituent  names. 


23  Blank  Card  to  end  the  flow  field 


USE  CARD  24  ONLY  IF  : 

ISCAL'A  = 

1 

24  12E6.3  1-6 

ABSCLI 

P> 

7-12 

ABSCLI 

!z) 

13-18 

ABSCLI 

I3 

19-24 

ABSCLI 

25-30 

ABSCLI 

is) 

31-36 

ABSCLI 

I6 

37-42 

DINSCI 

1 

43-48 

DINCCI 

2 

49-54 

DINSC(3) 

55-60 

DINSCI 

4) 

61-66 

DINSCI 

5) 

67-72 

D INSC (6 ) 

USE  CARDS  25  THRU  31 

ONLY  IF 

IO 

(Absorption  coeff.  scale  factors. 
[The  number  order  is  the  same  as  the 
lorder  in  the  IDENTC  list. 


Fine  structure  parameter  scale  factors 


= 1 (Absorption  coefficient  to  be  modified) 


25 


8F10.6 


26 


715 


1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

5 


G1 

62 

63 

G4 

G5 

G6 

IM0D 


Broadening  parameter  to  be  used  for  un- 
identified gases  in  the  flow  field. 
Number  order  is  the  same  as  the  order  in 
the  IDENTC  list.  For  example,  the  un- 
identified foreign  gas  broadening 
parameter  for  C02  is  G2. 

Number  of  constituents  for  which  absorp- 
tion coefficients  will  be  modified.  (If 
only  6 values  on  Card  25  are  to  be 
changed,  this  should  be  zeros). 


10 

NUMBER (1 ) 

NUMBER  is  the  position  of  the  constituent 

15 

A 

i 

to  be  modified  in  the  IDENTC  list.  For 

20 

25 

30 

y 

r 

example,  for  C02  NUMBER  (IM0D)  = 2. 

35 

NUMBER(IMOD) 
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INPUT  CARDS  (Continued) 


Card  Format  Col 

Computer 

Name 

Units 

Description 

Use  IM0D  sets  of  cards 

27  thru  31. 

The  sets  should  be  put  in  the  same  sequence 

as  used  in  the  NUMBER 

sequence  on 

card  26. 

27  815  1-5 

NUL(I) 

NUU(I) 

Initial  wavenumber  In  the  table. 

6-10 

Final  wavenumber  in  the  table. 

11-15 

I DNU  ( I ) 

Wavenumber  increment  In  the  table. 

20 

ICO 

0 absorption  coefficients  will  not 
be  modified. 

1 absorption  coefficients  will  be 
modified. 

25 

IFI 

0 fine  structure  parameters  will  not  be 
modified. 

1 fine  structure  parameters  will  be 
modified. 

USE  CARDS  28  AND  29  ONLY  IF  IC0>0 


28 


15, 5X 
7E10.6 


5 

11-20 

21-30 

31-40 

41-50 

51-60 

61-70 

71-80 


NACT(I) 


TEMAC ( I , K) 


Number  of  temperatures  In  table  (7  max) 
Temperatures  in  the  absorption  coefficient 
table  in  order  from  lowest  to  highest. 

K = 1,  NACT(I). 


Number  of  cards  29  is  (NUU(I)  - NUL ( I ) )/ I DNU ( I ) + 1 


29 


8E10.6 


1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

61-70 

71-80 


ENUK(L)  cm"1  Wavenumber  for  these  absorption  coefficients 
C0EF(k,L)  (cm-1  Ltd  Absorption  coefficients  for  each  temperature 
▲ in  the  table  (K  = 1,  NACT(I)). 


* 

C0EF(K,L) 


Use  cards  30  and  31  only  if  IFI  >0 


30 


15, 5X 
7E10.6 


5 NFINT(I) 

11-20  TEMFIN(I,K)  K 
21-30  ▲ 

31-40 

41-50 

51-60 

61-70  ▼ 

71-80  TEMFIN(I.K) 


Number  of  temperatures  in  table  (7  max). 
Temperatures  in  the  table  of  line 
densities  in  order  from  lowest  to 


highest,  K = 1 , NFINT(I) 
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INPUT  CARDS  (Concluded) 


Card 

Format 

Col 

Computer 

Name 

Units 

Description 

Number  of  cards 

31  is 

(NUU(I)  - NUL ( I ) )/ I DNU ( 1 ) +1 

31 

8E10.6 

1-10 

ENUK(L) 

„ -i 

cm 

Wavenumber  for  these  line 

11-20 

COEF (1 ,L) 

cm 

parameters.  Line  density 

21-30 

▲ 

each  of  the  temperatures. 

31-40 

NFINT(I) 

41-50 

51-60 

61-70 

1 

r 

71-80 

COEF(K.L) 

USE  CARD  32  ONLY  IF  I PLOT  >0 


32 


1615 


3-5 

8-10 


4 


JPLOT(I) 
1=1 ,NPLOT 


78-80 

3i5 

etc. 


Digits 


10 


20 


A series  of  NPLOT  three  digit  numbers 
will  be  used  to  give  plotting  instructions. 
The  first  digit  will  specify  the  units  of 
the  abscissa:  1 = wavenumber  - 1/cm; 

2 = wavelength  - microns.  The  remaining 
two  digits  will  specify  the  variable  on 
the  ordinate  as  follows: 


Variable 

Name 


Description 


Spectral  irradiance  in  watts/cm2-y-sr 

(NSPEK=1  or  2) 

" "watts/cm2-y  (NSPEK=3  or  4) 

" "watts/cm2 -sr-cm_1  (NSPEK  = 11  or  12 
" "watts/cm2-cnf 1 (NSPEK=13  or  14) 

Transmissivity 


31  thru  37  FO 


41  thru  46  ACO 


51  thru  56  ADO 


F # 


Optical  depth  for  the  weak  line  limit, X*. 

Last  digit  indicates  constituent  according 
to  IDENTC  list. 

Collision  broadened  fine  structure  parameter, 
parameter,  a^  Last  digit  indicates 

constituent. 

Doppler  brwdended  fine  structure 
parameter  aD.  Last  digit  indicates 

constituent. 
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INPUT/OUTPUT  TAPES 

The  use  of  tapes  by  the  program  depends  upon  the  options  selected.  It 
can  operate  without  input  or  output  tapes,  or  it  can  prepare  output  tapes 
which  may  later  be  used  as  input.  Each  tape  will  be  described  as  its 
primary  function  with  notes  to  explain  the  optional  use. 

Flow  Field  Tape 

This  tape  is  used  when  KASE  is  positive  and  I FLOW  is  2 or  3.  It  pro- 
vides input  gas  property  data.  The  Lockneed  developed  programs  (Refs.  1 
and  2)  can  prepare  an  axi symmetric  flow  field  tape,  but  no  method  is  presently 
available  for  automated  preparation  of  a three-dimensional  flow  field  tape. 

Any  flow  field  input  on  cards  or  tape  can  be  transferred  to  a tape  using 
the  IFLOTP  option. 

Density  - 800  BPI  FORTRAN  IV  Binary 
For  IFLOTP  = 0:  input  on  Logical  8 
For  IFL0TP>0:  input  on  Logical  11  (If  KASE  >0) 

output  on  Logical  8 

The  composition  of  the  input  tape  varies  slightly  depending  upon  whether 
the  flow  field  is  axisymmetric  or  three-dimensional.  Therefore,  the  tapes 
will  be  described  separately.  Several  flow  field  cases  may  be  stacked  on 
one  tape.  The  last  case  on  a tape  is  followed  by  an  end  of  file, or  a record 

in  which  each  of  the  twelve  words  in  HDG(I)  is  ******. 

For  I FLOW  = 2 

For  the  axisymmetric  flow  field,  the  first  three  records  consist  of 

title  information  and  constituent  identification: 

Record  1 Words  1,  12  (HDG(I), 1=1 ,12)  Flow  field  problem 

title.  For  information 
only. 
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INPUT/OUTPUT  TAPES  (Continued) 

Record  2 

Words  1 ,4 

(BETA(I) ,1=1 ,4 

Propellant  or  gas 
title.  (For  infor- 
mation only.) 

Record  3 

Word  1 

N0C0N 

Number  of  constituents 
on  the  flow  field  tape. 

Words  2,3 

(KONAME(I.J), 
0=1 ,2) 

Name  of  the  Ith 
constituent.  Consists 
of  two  A6  words. 

Word  4 

CONWT(I) 

Constituent  molecular 
weight. 

Word  5 and 

subsequent  - Repeat  words  2 thru  4 for  I = 1 , N0C0N 

The  remainder  of  the  case  consists  of  two-record 
gas  data. 

groups  specifying  the 

Record  4 

Word  1,3 

NOPS 

Number  of  radial  values 
for  property  data  at 
this  Z-plane.  This  is 
repeated  3 times. 

Word  4 

LS(I) 

When  ISHOCK  =-l , this 
value  must  be  the  (L) 
index  of  the  radial 
point  just  inside  the 
shock.  If  ISHOCK  f -1, 
N0PS  may  be  repeated  here 
to  maintain  proper  record 
length. 

Word  5 

ZZ(I) 

Z-coordinate  of  the  data 
presented  in  the  next 
record. 

Record  5 

Word  1 
Word  2 
Word  3 
Word  4 

R(I,L) 
T(I.L) 
P(I.L) 
FF(L,J), 
0=1 ,N0C0N 

Radius 

Temperature 

Pressure 

Mole  fraction.  Repeated 
for  J=1 ,N0C0N. 

The  word  group  above  (Record  5)  is  repeated  on  this  record 

for  L=1,N0PS. 

Sets  of  records  like  4 and  5 are  repeated  until  all  the  data  is  loaded.  The 
end  of  the  flow  field  is  signaled  by  a record  with  NOPS  = 0. 
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INPUT/OUTPUT  TAPES  (Continued) 


For  I FLOW  * 3 

For  the  three-dimensional  flow  field,  the  first  three  records  consist  of 


title  information,  flow  field  geometric  parameters,  and  constituent  identi- 
fication. 


Record  1 

Words  1,12 

(HDG(I),I  = 1,12) 

Flow  field  problem 
title.  For  infor- 
mation only. 

Record  2 

Words  1 , 4 

(BETA(I).I  = 1,4) 

Propellant  or  gas  title. 
(For  information  only). 

Word  5 

ISYMSC 

Flow  field  geometry 
option. 

Word  6 

HANG 

Half  angle  of  symmetri- 
cal sector  used  in  flow 
field  geometry. 

Word  7 

CTCD 

Outboard  engine  spacing 
from  vehicle  centerline. 

Record  3 

Word  1 

N0C0N 

Number  of  constituents 
on  the  flow  field  tape. 

Words  2,3 

(K0NAME(I,J),J=1 ,2)Name  of  the  Ith  consti- 
tuent. Consists  of  two 
A6  words. 

Word  4 

CONWT(I) 

Constituent  molecular 
weight. 

Word  5 and  subsequent  - Repeat  words  2 thru  4 for  I = 1 ,N0C0N 

The  remainder  of  the  case 
the  gas  data. 

consists  of  groups  of  three  records  specifying 

Record  4 

Word  1 

zz(i) 

2 coordinate  of  the  data 
in  the  following  records. 

Word  2 

NETAS 

Number  of  n -values  to 
be  specified  at  this  Z- 
plane. 
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Record  5 and  6 are  repeated  for  0=1,  NETAS 


Record  5 

Word  1 

ETA(I,J) 

Word  2 

NOPS 

Word  3 

LS(I,J) 

Record  6 

Word  1 
Word  2 
Word  3 
Word  4 

R(U,L) 
T I.J.L 
P(I,0,L) 
FF(L,N),N  = 1 

RTR  013-1 


Value  of  n at.  the  Jth 
position. 

Number  of  radial  points 
for  gas  property  values 
at  this  n -position. 

When  ISHOCK  = -1  this 
value  must  be  the  (L) 
index  of  the  radial 
point  just  inside  the 
shock.  If  ISHOCK /-l , 
NOPS  may  be  repeated  here 
to  maintain  proper  record 
length. 

Radi  us 
Temperature 
Pressure 
Mole  fraction. 

Repeated  for  N = 1 , 

N0C0N. 


The  word  group  (Record  6)  above  is  repeated  on  this  record 
for  L = 1,  NOPS. 


Sets  of  records  like  4 thru  6 are  repeated  until  all  the  data  is  loaded. 
The  end  of  the  flow  field  is  signaled  by  a record  with  NETAS  = 0. 


Line  of  Sight  Tape 

Density  - 800  BPI  FORTRAN  IV  Binary  Input  or  output  on  Logical  10 
This  tape  is  generated  by  the  flow  field  subroutines  and  can  be  used  as  an 
input  tape  with  IFLOW  = 1.  The  primary  purpose  for  saving  this  tape  is  to 
allow  restart  and  completion  of  a radiation  calculation  that  has  been 
terminated  by  the  estimated  run  time  limit.  It  can  also  be  used  to  rerun 
a problem  with  reduced  integration  limits.  However,  in  this  case  DS  must  not 
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INPUT/OUTPUT  TAPES  (Continued) 


be  changed  and  the  values  of  0 and  $ for  the  lines  of  sight  must  match 
those  existing  on  the  tape.  The  first  several  records  consist  of  input 
data  which  is  defined  in  previous  Section  on  input  cards. 

Logical  record  one  consists  of  37  + K'  CON  words  in  the  following 
order: 

Word  1 KASE  Word  15  SMAX  Word  29  TANGLE 


Word  29  TANGLE 


2 DATE 

16  DS 

30  ENUI 

3 I FLOW 

17  XP 

31  ENUF 

4 HDG 

18  YP 

32  DENU 

5 BETA 

19  ZP 

33  TDIFF 

6 IUNITS 

20  I PLANE 

34  ACMIN 

7 IZI 

21  NOKON 

35  SZERO 

8 IZF 

22  IBLOCK 

36  (NUMCON(K), 

9 THETAI 

23  ICARB 

K = 1 , NOKON) 

10  THETAF 

24  I SHOCK 

11  DIHETA 

25  ISCALF 

12  PHII 

26  HANG 

Word  (36  + NOKON)  ISYMSC 

13  PHIF 

27  CTCD 

Word  (37  + NOKON)  KCARB 

14  DPHI 

28  RSCAL 

The  geometric  scaling  data 

is  listed  in  the  second  record.  It  consists  of  9 

words  as  follows: 

Word  1 RSCAL 

Word  4 FISCAL 

Word  7 F4SCAL 

2 TSCAL 

5 F2SCAL 

8 F5SCAL 

3 PSCAL 

6 F3SCAL 

9 F6SCAL 
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INPUT/OUTPUT  TAPES  (Continued) 

Logical  record  three  gives  the  orientation  of  the  plane  of  interest  and  will 
consist  of  either  one  or  nine  words  depending  upon  IPLANE. 

If  IPLANE  = 0 
Word  1 OMEGA 
If  IPLANE  = 1 


Word  1 ULX 

Word  4 VLX 

Word  7 

WLX 

2 ULY 

5 VLY 

8 

WLY 

3 ULZ 

6 VLZ 

9 

WLZ 

Logical  record  four  lists 

the  constituent  names 

used  for 

this 

constituent  name  is  two  A6  words. 

Word  1 

KON  (1) 

(2*N0K0N) 

KON  (2*N0K0N) 

The  remainder  of  the  input  data  may  require  from  1 to  173  additional  records 
depending  upon  the  input  options  chosen.  Up  to  151  records  may  be  required  for 
specifying  the  addition  of  carbon  particles  to  the  plume  (ICARB>0).  The 
blocking  circle  specifications  requires  one  record  if  IBL0CK;  0,  and  the 
specification  of  shock  radii,  for  ISHOCK  =1,  may  require  up  to  21  records. 

The  arrangement  of  the  records  and  the  dependence  upon  the  options  chosen 
are  outlined  below. 

The  number  of  records  required  for  adding  carbon  particles  to  the  flow 
field  depends  upon  the  values  of  ICARB  and  I FLOW. 

If  ICARB  = 1,  one  record  with  two  words  is  used. 

Word  1 FC  (1,1,1) 

2 CRATIO  (1,1,1)  - This  word  not  required,  but  is  written  for 

convenience 
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INPUT/OUTPUT  (Continued) 

If  ICARB  = 2,  two  to  eight  records  may  be  used. 


Record  1 

Word  1 

ZZFC  (1) 

This  is  not  required  under  this  option, 
but  is  written  here  to  simplify 
routines. 

2 

NOETAC  (1) 

For  IFLOW  = 2,  this  is  set  equal 
to  1 

Record  2 

Word  1 

ETAC  (1 ,L) 

This  is  required  for  IFLOW  = 3 only 

2 

NOPSFC  (1,L) 

Record  3 

Word  1 

FC  (1 ,L,K) 

Sets  of  FC  and  CRATIO  are  repeated 
for  K = 1,  NOPSFC  (1,0 

2 

CRATIO  (1 ,L,K) 

Records  2 

and  3 are  repeated  for  L - 

1,  NOETAC (1). 

ICARB  = 3,  i 

up  to  151 

records  may  be  used. 

Record  1 

Word  1 

NZFC 

Record  2 

Word  1 

ZZFC ( I ) 

2 

NOETAC ( I ) 

Set  = 1 for  IFLOW  = 2 

Record  3 

Word  1 

ETAC ( I ,L ) 

Required  for  IFLOW  = 3 only,  but 
written  in  all  cases  for  simplicity 

2 

NOPSFC  (1,0 

Record  4 

Word  1 

FC(I,L,K) 

Sets  of  FC  and  CRATIO  are  repeated 
for  K = 1,  NOPSFC ( I ,L) 

2 

CRATIO(I,L,K) 

Records  3 

and  4 are  repeated  for  L = 

1,  NOETAC  (I) 

Records  2 

thru  4 are  repeated  for  I ! 

= 1 , NZFC 
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The  record  for  the  blocking  circle  specifications  consists  of  IBLOCK 
groups  of  6 words  each. 

Word  1 XC(I) 

2 YC  ( I ) 

3 ZC(I) 

4 RC(I) 

5 KIND ( I ) 

This  word  group  is  repeated  for  1=1,  IBLOCK. 

The  number  of  records  required  for  listing  the  shock  radii  when  ISHOCK  = 1 
depends  upon  the  flow  field  geometry.  If  IFLOW  = 2,  only  two  records  are 
used,  but  if  IFLOW  = 3,  up  to  21  records  may  be  required. 

Record  1 Word  1 NZSH 

2 IZSH 

3 N0ETAS(I ) This  word  is  repeated  for  I = 1 , 

NZSH.  This  information  is  not 
required  for  IFLOW  = 2,  but  words 
are  placed  on  the  record  as  shown. 

If  IFLOW  = 2 only  one  more  record  is  required 

Record  2 Word  1 RSHOCK  (1,1)  This  word  is  repeated  or  I = 1 , 

NZSH 

If  IFLOW  = 3 the  number  of  additional  records  is  equal  to  NZSH. 

Record  2 Word  1 RSH0CK(I,J)  This  word  is  repeated  for  J 3 1, 

NOETAS(I) 

Record  2 is  repeated  for  1=1,  NZSH. 

This  concludes  the  records  specifying  input  data.  Following  these  records, 
are  a number  of  records  equal  to  the  number  of  lines  of  sight  in  the  problem 
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INPUT/OUTPUT  TAPES  (Concluded) 

plus  1.  The  final  record  signals  the  end  of  the  problem  by  having  LSET  = 0. 
The  words  in  each  record  are  as  follows: 


Word  1 

LSET 

Number  order  of  the  line  of  sight 

2 

LPT 

Number  of  points  on  the  line  of  sight 

3 

L BLOCK 

Blocking  index 

4 

SINKO 

Product  of  sine  cose 

5 

THETA 

Angle  6 

6 

PHI 

Angle  <f> 

7 

DSFACT 

Differential  angle  factor  (sine  cose  ded<j>/Tr) 

8 

DX 

Integration  increment  on  X coordinate 

9 

DY 

Integration  increment  on  Y coordinate 

10 

DZ 

Integration  increment  on  Z coordinate 

11 

SMIN 

Upper  limit  of  the  s summation 

12 

XT  (L ) 

X coordinate 

13 

YT(L) 

Y coordinate 

14 

ZT(L) 

Z coordinate 

15 

TWT(L) 

Temperature 

16 

PWT(L) 

Pressure 

17  to  (16+N0K0N)FWT(L,K)  Mole  fraction.  This  is  repeated  K = 1,  N0K0N 

Words  12  thru  (16+NOKON) 

are  repeated  for  L = 1 , LPT. 

Final  Word 

DPHI 

Integration  increment  for  <j> 
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PRINTED  OUTPUT 

The  initial  section  of  output  is  printed  in  the  FLOUT  subroutine.  It 
consists  of:  the  problem  and  flow  field  identification,  the  unit  system  for 
input,  the  limits  of  integration,  and  the  coordinates  of  the  point  of 
interest.  Following  this,  the  inclination  angle  OMEGA  or  the  nine  directional 

cosines  (ULX  WLZ)  are  output  to  define  the  orientation  of  the  plane 

containing  the  point  of  interest.  Then  the  constituents  which  will  be  con- 
sidered are  listed.  After  this,  the  output  depends  upon  the  options  chosen. 

If  ICAR8  is  greater  than  0,  the  desired  carbon  particle  distribution  is  listed. 
Following  this,  the  blocking  circle  specifications  are  listed  if  IBLOCK  Is 
greater  than  0.  Then  the  shock  radii  data  are  output  if  ISHOCK  is  greater 
than  0.  Finally,  this  section  of  output  is  completed  by  printing  the 
geometric  scale  factors  if  ISCALF  is  greater  than  0. 

The  second  section  of  output  Is  printed  in  FL0W2D  or  FL0W3D  if  IFLOW  is 
greater  than  1.  If  IFLOW  is  less  than  or  equal  to  1,  this  section  of  output 
is  omitted.  The  output  quantity  depends  upon  the  option  NFL0W.  If  NFLOW  = 1, 
the  radius  (R),  temperature  (T),  pressure  (P),  and  constituent  mole  fractions 
(FI  to  F6)  at  each  point  In  the  flow  field  are  listed.  If  NFLOW  = 0,  only  the 
first  and  last  radial  points  are  listed  to  assure  proper  input  operation. 

The  third  section  of  output  is  printed  in  ACDATA.  This  section  lists  the 
absorbtion  coefficient  and  line  density  data  which  will  be  used  in  the 
program.  If  NABS  = 0,  only  the  first  and  last  entry  for  each  constituent 
are  listed  to  minimize  output  while  assuring  proper  program  operation  . If 
NABS  = 1,  all  values  for  each  constituent  will  be  listed. 
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PRINTED  OUTPUT  (Continued) 

The  fourth  section  of  output  Is  printed  in  either  M0DL3A  or  M0DEL3, 
depending  on  the  value  of  IRAD.  The  output  depends  upon  the  option  NLINE. 
If  NLINE  - 0,  the  integrated  irradiance  (q/A  or  FLUX)  is  printed  after  the 
integration  over  each  A9  step  is  complete.  When  NLINE  = 1 , the  output  is 
printed  for  each  line  of  sight.  In  this  case,  the  value  of  6 (THETA)  and 
4)  (PHI)  for  each  line  of  sight  is  listed  along  with  the  following  data: 


OUTPin  LABEL  MNEMONIC* 

SHAPE  FACTOR  - DELTA  DSFACT 

SUM  SFACT 

FLUX  - DELTA  LOSFLX 

SUM  FLUX 

RADIANCE  - LINE  OF  SIGHT  RADLOS 

- AVERAGE  RADAVG 


In  addition  to  the  above  data,  the  output  will  note  which  lines  of  sight 
missed  the  gas  plume  and  which  lines  of  sight  were  blocked  by  blocking  circles 
before  entering  the  gas  plume.  For  the  lines  of  sight  which  miss  the  gas,  DX,  DY, 
and  DZ  are  printed  to  Indicate  the  direction  of  the  line  of  sight. 

The  output  in  the  fourth  section  is  similar  for  NLINE  = 2 or  3.  In  both 
cases  the  coordinates  for  the  position  on  the  line  of  sight  are  listed  along 
with  the  gas  properties,  number  of  As  steps,  incremental  and  summed  irradiance, 
and  average  transmissivity.  The  difference  between  NLINE  = 2 or  3 is  that  for 
NLINE  = 2 the  output  is  printed  at  the  end  of  each  Integration  step  (DSS), 


* MNEMONICS  are  defined  in  APPENDIX  A. 
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PRINTED  OUTPUT  (Continued) 

while  with  NLINE  = 3,  the  output  is  printed  at  each  geometric  step  (DS)  along 
with  the  average  values  for  each  Integration  step.  Due  to  the  storage  allocation 
system  used  in  M0DE1.3,  only  the  data  for  the  first  spectral  interval  along  the 
line  of  sight  will  be  printed.  The  computer  names  corresponding  to  the  output 
labels  are  listed  below: 


OUTPUT  LABEL 
X 

Y 

Z 

X 

T 

P 

FI  thru  F6 
DS  STEPS 
FLUX-FOR  DS 
-SUM 
AVG  - G 


MNEMONIC 

X 

Y 

Z 

X 

TW 

PW 

FW(1 ) thru  FW(6) 

JFK 

DSFLUX 

LOSFLX 

AVGG 


At  the  beginning  of  the  line  of  sight,  the  output  values  of  6 and  <{>  for 
the  line  of  sight  are  listed,  and  at  the  end  of  the  line  of  sight  the  following 


data  is  printed. 

OUTPUT  LABEL  MNEMONIC 

SHAPE  FACTOR  - PLUME  SFACT 

TOTAL  TSFACT 

FLUX  - TOTAL  FLUX 

LOS  LOSFLX 

RADIANCE  - AVG  RADAVG 


78 


REMTECH  INCORPORATED 


RTR  013-1 


OUTPUT  LABEL  MNEMONIC 

PATH  LENGTH  FOR  H20  UH20 

C02  UC02 

CO  UCO 

Cl  UCARB 

HCL  UHCL 

HF  UHF 


When  the  integration  is  complete,  the  following  output  is  printed 
regardless  of  the  value  of  NLINE: 


OUTPUT  LABEL 

MNEMONIC 

NUMERICAL  PLUME  SHAPE  FACTOR 

SFACT 

NUMERICAL  TOTAL  FACTOR 

TSFACT 

EXACT  TOTAL  SHAPE  FACTOR 

TRUESF 

PLUME  RADIANCE 

RADPLM 

OVERALL  RADIANCE 

RADAVG 

TOTAL  FLUX 

FLUX 

The  final  output  item  in  the  fourth  section  is  the  calculation  time  for 
the  radiation  calculation. 

The  fifth  section  of  output  data  is  printed  in  SPCOUT.  It  provides  for 
spectral  output  and  plotting.  The  printed  spectral  data  is  output  when 
NSPEK  >0.  In  the  spectral  output,  some  of  the  output  data  is  modified  to 
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PRINTED  OUTPUT  (Continued) 

different  units.  Therefore,  some  of  the  output  items  listed  below  will  be 
redefined. 


OUTPUT  LABEL 

MNEMONIC 

DEFINITION 

WAVENUMBER 

ENU 

Wavenumber  - 1/cm 

WAVELENGTH 

ELAM 

Wave! ength-mi crons 

HoVenu) 

IRRADIANCE-WATTS/CM-STER 

SFLXTO(L) 

Spectral  radiance  per 
unit  wavenumber.  NSPEK  = 
11  or  12 

-WATTS/CM 

SFLXTO(L) 

Spectral  irradiance 
per  unit  wavenumber 
NSPEK  = 13  or  14 

-WATTS/SQ-CM-MICRON-STER 

SFLXTO(L) 

Spectral  radiance  per 
unit  wavelength.  NSPEK  = 
1 or  2 

-WATTS/SQ-CM-MICRON 

SFLXTO(L) 

Spectral  irradiance 
per  unit  wavelength. 
NSPEK  = 3 or  4 

The  following  data  is  for  the  last  line 

of  sight  only. 

G 

GNEW(L) 

Transmissivity 

OPTICAL  DEPTH 

0D(I) 

Optical  depth  for  all 
constituents  -math 
symbol  X(s,w). 

The  following  data  is  output  only  if  NSPEK  = X2  or  X4. 

CONSTITUENT  IDENTIFICATION 

KON(I) 

Constituent  name  for 
data  on  this  line 

DSTAR 

F0(I)  or 
DSTAR 

Optical  depth-math 
symbol  X*. 

D/DSTAR 

DODST 

Ootical  depth  ratio  for 
the  constituent  indicated 
•math  symbol  X/X* 
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OUTPUT  LABEL 

MNEMONIC 

DEFINITION 

DC/DSTAR 

DCODST 

Collision  broadened 
optical  depth  -math 
symbol  Xc/X* 

DD/DSTAR 

DDODST 

Doppler  broadened  optical 
depth  -math  symbol  XD/X* 

AC 

AC 

Collision  broadened  fine 
structure  parameter  -math 
symbol  ac 

AD 

AD 

Doppler  broadened  fine 
structure  parameter  -math 
symbol  a 

When  IRAD>0  and  NSPEK  = X2  or  X4,  the  detailed  spectral  output  must  be  limited 
because  of  storage  limitations.  This  will  occur  when  (ENUF-ENUI)/DENU  * IRAD 
> NUDIM  (NUDIM  is  set  at  840).  If  the  storage  limit  is  exceeded,  only  a 
portion  of  the  spectral  data  at  the  larger  values  of  ENU  will  be  output. 


i 

i 


f 
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DIAGNOSTICS 

If  the  number  of  line  groups  (IRAD)  exceeds  10,  the  program  will  stop  and 
print 

NUMBER  OF  LINE  GROUPS  EXCEEDS  THE  MAXIMUM  ALLOWED  VALUE  OF  10. 

If  the  number  of  spectral  intervals  exceeds  440,  the  program  will  stop  and 
print 

(ENUF-ENUI)/IDENU+1  IS  GREATER  THAN  THE  MAXIMUM  NUMBER  OF  SPECTRAL  INTERVALS 
DIMENSIONED  (440)  . CASE  TERMINATED. 

If  the  flow  field  case  to  be  read  from  tape  cannot  be  found  on  the  tape, 
the  program  will  return  to  read  the  next  data  and  write 

*****S0RRY,  CASE  XXXXXXXXXX  IS  NOT  ON  THIS  TAPE.  RUN  TERMINATED. 

If  the  flow  field  from  tape  does  not  contain  the  desired  constituents,  the 
program  returns  to  read  new  data  and  writes 

*****CASE  NUMBER  XXXXXXXXXX  DOES  NOT  CONTAIN  THE  CONSTITUENT  'XXXXXXXXXXXX' 

If  an  input  constituent  cannot  be  identified  by  the  program,  the  program 
returns  to  read  new  data  and  writes 

THE  CONST ITUENT-XXXXXXXXXXXX- IS  NOT  RECOGNIZED  BY  THE  PROGRAM.  CHECK 
AGAINST  NAMES  IN  THE  IDENTC  LIST. 

If  NZFC  exceeds  10,  the  program  will  return  to  read  new  data  and  write 
NUMBER  OF  CARBON  Z-CUTS  EXCEEDS  10 

If  NOETAC(I)  exceeds  7,  the  program  will  return  to  read  new  data  and  write 

NUMBER  OF  CARBON  ETA-CUTS  EXCEEDS  7.  NOETAC i VXX )=XXX 

If  N0PSFC(I,L)  exceeds  10,  the  program  will  return  to  read  new  data  and 
write 

NUMBER  OF  CARBON  RADIUS  POINTS  EXCEEDS  10l  N0PSFC(XXX,XXX)=XXX 

If  the  number  of  Z-planes  for  shock  radii  exceeds  20,  the  program  will 
return  to  read  new  data  and  write 

NUMBER  OF  SHOCK  RADIUS  Z-CUTS  EXCEEDS  20. 
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If  the  number  of  n-values  for  shock  radii  exceeds  13,  the  program  will 
return  to  read  new  data  and  write 

NUMBER  OF  SHOCK  RADIUS  ETA-CUTS  EXCEEDS  13.  N0ETAS(XXX)=XXX 

If  IFL0W=1  and  the  input  value  of  SMAX  exceeds  that  on  tape  10,  the  program 
will  return  to  read  new  data  and  write 

INPUT  SMAX  XXXXX.XX  EXCEEDS  LOS  TAPE  VALUE  OF  XXXXX.XX 

If  IFL0W=1  and  the  input  value  of  DTHETA  does  not  agree  with  that  of  tape 
10,  the  program  will  return  to  read  new  data  and  write 

INPUT  DTHETA  XXXXX.XX  NOT  EQL  LOS  TAPE  VALUE  OF  XXXXX.XX 

If  IFL0W=1  and  the  input  value  of  D?HI  does  not  agree  with  that  on  tape 
10,  the  program  will  return  to  read  new  data  and  write 

INPUT  DPHI  XXXXX.XX  NOT  EQL  LOS  TAPE  VALUE  OF  XXXXX.XX 

If  IFL0W=1  and  the  0 limits  input  exceed  those  from  tape  10,  the  program 
will  return  to  read  new  data  and  write 

THETA  LIMITS  DO  NOT  AGREE  - INPUT  VALUES  ARE  THETAI=XXXXX.XX  AND  THETAF= 
XXXXX.XX  LOS  TAPE  VALUES  ARE  THETAI=XXXXX.XX  AND  THETAF=XXXXX.XX 

If  IFL0W=1  and  the  <t>  limits  input  exceed  those  from  tape  10,  the  program 
will  return  to  read  new  data  and  write 

PHI  LIMITS  DO  NOT  AGREE  - INPUT  VALUES  ARE  PHI I=XXXXX . XX  AND  PHIF=XXXXX.XX 
LOS  TAPE  VALUES  ARE  PHI I=XXXXX . XX  AND  PHI F=XXXXX . XX 

If  the  number  of  line  of  sight  points  exceeds  the  maximum  allowed,  the 
program  will  return  to  read  new  data  and  write 

NO.  OF  LINE  OF  SIGHT  POINTS  EXCEEDS  XXXX,  NEED  LARGER  DELTA  S 

If  it  is  necessary  to  delete  radial  flow  field  properties  to  fit  within 
storage  limitations,  the  program  will  write  the  following  warning  and 
continue 

*****RADIAL  POINTS  WERE  DELETED  FROM  THE  FOLLOWING  DATA.  NRMAX=XXXXX 
NOPS=XXXXX 

If  the  flow  field  exceeds  the  maximum  Z-plane  at  which  the  carbon  mole 
fraction  is  specified  when  ICARB=3,  the  carbon  mole  fraction  will  be  omitted, 
and  the  program  will  write 

**********N0  CARBON  MOLE  FRACTION  SPECIFIED  FOR  Z=XXXXX.XX 
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If  the  number  of  n-values  exceeds  13  when  IFL0W=3,  the  program  will  stop 
and  write 

NUMBER  OF  ETA-CUTS  EXCEEDS  13.  NETA(XXX )=XXX 

If  the  absorption  coefficients  are  being  modified  and  the  number  of  entries 
in  the  coefficient  table  (COEF ( I s K) ) exceeds  the  maximum  allowed,  the  pro- 
gram returns  to  read  the  next  case  and  writes 

LAST  COEF  INDEX  OF  XXXXX  EXCEEDS  THE  MAXIMUM  ALLOWED  VALUE  OF  XXXXX 

If  the  optical  depth  is  negative,  the  transmissivity  is  taken  to  be  equal 
to  that  at  the  previous  line  of  sight  increment  and  the  program  writes 
the  following  warning 

*****NEGATIVE  DM0  AT  NU=XXXXXX.  DMO=XX.XXXXE+XX 

If  the  temperature  is  negative  at  any  point,  the  following  warning  will  be 
wri tten 

******WARNING  . NEGATIVE  TEMPERATURE  ENCOUNTERED  AT  LSET=XXXXX  LOSPT=XXXXX 

If  the  number  of  spectral  intervals  summed  over  all  the  constituents  exceeds 
NUDIM,  the  program  returns  to  read  new  data  and  writes 

*****THE  SUM  OF  SPECTRAL  INTERVALS  FOR  ALL  CONSTITUENTS  IS  XXXX.  THE  MAX- 
IMUM ALLOWED  IS  XXXX. 

In  order  to  assure  that  the  plot  section  of  the  program  was  reached,  the 
following  statements  will  be  printed  at  appropriate  locations 

BEGINNING  PLOT  SECTION 
END  OF  PLOT  SECTION 

If  the  abscissa  indication  is  missing  for  a plot,  the  program  will  skip 
to  the  next  plot  and  write 

ABSCISSA  INDICATION  MISSING 

If  the  ordinate  indication  is  missing  for  a plot,  the  program  will  skip 
to  the  next  plot  and  write 

ORDINATE  INDICATION  MISSING 

If  an  input  error  calls  for  a plot  of  the  fine  structure  parameters  for 
carbon  particles,  the  plot  will  be  skipped  and  the  following  message  will 
be  written 

FINE  STRUCTURE  PARAMETERS  NOT  USED  FOR  CARBON.  PLOT  XXX  DELETED 
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PLOTS  GENERATED 

The  number  of  plots  generated  will  be  equal  to  the  Input  option 
NPLOT.  The  data  plotted  and  the  ordinate  and  abscissa  descriptions 
will  vary  depending  upon  the  specifications  set  forth  by  NSPEK  (input 
card  2)  and  JPLOT  (input  card  32).  The  first  digit  of  JPLOT  specifies 
the  abscissa  variable  and  the  second  two  digits  specify  the  ordinate 
variables.  The  arrangement  used  is  defined  in  the  table  below: 


Abscissa  Selection  (First  Digit) 


JPLOT  Abscissa 

1XX  Wavenumber  (ENUV(I)) 

2XX  Wavelength  (ELAMV(I)) 


Ordinate  Selection  (Second  and  Third  Digits) 


JPLOT 
XI 0 


X20 

X3X 

X4X 

X5X 


NSPEK  Ordinate 

01  or  02  Spectral  radiance  per  unit  wavelength 

watts/cm2-y-sr 

03  or  04  Spectral  irradiance  per  unit  wavelength 

watts/cm2-y 

11  or  12  Spectral  radiance  per  unit  wavenumber 

watts/cm2-sr-cm“ 1 

13  or  14  Spectral  irradiance  per  unit  wavenumber 

watts/cm2-cm" 1 

Transmissivity  (GNEW(L)) 

Optical  depth  for  the  weak  line  limit 

mm 

Collision  broadened  fine  structure  parameter 
(ac)  averaged  over  the  path. 

Doppler  broadened  fine  structure  parameter 
(aD)  averaged  over  the  path. 
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